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ABSTRACT 

With the goal of identifying high redshift radio galaxies with FRI classification, here are pre- 
sented high resolution, wide-field radio observations, near infra-red and optical imaging and 
multi-object spectroscopy of two fields of the Leiden-Berkeley Deep Survey. These fields, 
Hercules. 1 and Lynx. 2, contain a complete sample of 81 radio sources with Si.4ghz > 0.5 
mjy within 0.6 square degrees. This sample will form the basis for a study of the population 
and cosmic evolution of high redshift, low power, FRI radio sources which will be presented 
in Paper II. Currently, the host galaxy identification fraction is 86% with 1 1 sources remain- 
ing unidentified at a level of r'^25.2 mag (Hercules; 4 sources) or r'>24.4 mag (Lynx; 7 
sources) or .K;>20 mag. Spectroscopic redshifts have been determined for 49% of the sample 
and photometric redshift estimates are presented for the remainder of the sample. 

Key words: galaxies: active - galaxies: evolution - galaxies: photometry - galaxies: distances 
and redshifts - radio continuum: galaxies 



1 INTRODUCTION 

Radio-loud active galaxies that display extended jet emission from 
their central cores can be divided into two main types; Fanaroff & 
Riley class I and II (FRI and FRII; Fanaroff & Riley 1974). The 
galaxies of FRI type are 'edge-darkened' with the majority of their 
emission confined to their central regions and jets that flare out 
close to the nucleus. On the other hand, the FRII galaxies are 'edge- 
brightened' meaning the bulk of their emission originates from the 
hotspots at the ends of their highly collimated jets. The FRII galax- 
ies are the more luminous of the two classes and typically have 
fl78MHz > 10 24 " 25 WHz-'sr" 1 but there is significant overlap 
at the break luminosity. The FRIs and FRIIs have also been sug- 
gested as the unbeamed parent populations of BL Lac objects and 
flat spectrum quasars respectively (Jackson & Wall, 1999). 

The differences between the two FR classes are not confined 
to the lobe morphology. For instance, Zirbel & Baum (1995) found 
that the FRIIs produce 10-50 times more emission line luminosity 
than the FRIs at a particular radio core power. Additionally, optical 
observations by Owen & Laing (1989) found that the host galaxies 
of FRIs tended to be larger and more luminous than those of FRIIs, 
though later work by Ledlow & Owen (1996) suggested that this 
result was caused by a combination of sample selection effect, and 
only observing a small range in radio power. 

It is not yet clear whether the observed morphological differ- 
ences between the two FR classes are the result of fundamental 
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differences in the properties of the central engine (e.g. lower accre- 
tion rates in FRIs leading to advection dominated accretion flow, 
or a slower FRI black hole spin) as advocated by e.g. Baum et. al. 
(1995), or differences in the interactions of the jets with their en- 
vironments as suggested by the work of Gopal-Krishna & Wiita 
(2000) and Gawronski et al. (2006). This intrinsic/extrinsic ques- 
tion is of vital importance for understanding the relationship be- 
tween these objects; if the intrinsic difference model is correct then 
the FRIs and FRIIs are two discrete classes of object, however if the 
evidence suggests the other model is correct, the underlying prop- 
erties of the classes would be the same. In the latter case the two 
classes may simply represent different stages in the evolution of a 
radio galaxy, i.e. it starts out as a powerful, high-luminosity FRII 
and as it ages its jets become less powerful and it becomes an FRI 
(e.g. Willott et al. 2001). 

One of the key ways in which the differences between the two 
classes can be investigated is through their evolution with cosmic 
epoch since, if the extrinsic model is correct, then FRIs and FRIIs 
of the same luminosity should evolve in similar ways. FRIIs are 
known to undergo strong cosmic evolution, with density enhance- 
ments for the most luminous of a factor of 100-1000 out to redshifts 
of 1-2, compared to a factor of ~10 for the less luminous, sources 
(Wall, 1980); the behaviour of the FRIs is less clear. Low-redshift 
studies initially suggested that they had a constant space density 
(e.g. Jackson & Wall 1999; Willott et al. 2001) and this appeared to 
be confirmed at higher redshift by Clewley & Jarvis (2004). How- 
ever, they selected the FRIs in their sample using a luminosity cut; 
this could lead to FRIs being missed (particularly the more lumi- 
nous FRIs which may evolve the most) since the FRI/FRII break 
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luminosity is not fixed, but is a function of host galaxy magnitude 
(Ledlow & Owen 1996). It is clear therefore, that to define a ro- 
bust sample of distant FRIs, in order to obtain an accurate picture 
of their high redshift behaviour, radio morphological classification 
is a necessity. 

Determining the cosmic evolution of FRI radio galaxies is also 
important because of the impact that they may have on galaxy for- 
mation and evolution. Models of galaxy formation are increasingly 
turning to these objects to solve the problem of massive galaxy 
over-growth (e.g. Bower et al., 2006). It is predominantly the lower 
luminosity sources that provide the necessary feedback for this, 
(Best et al., 2006), and may possibly be limited to the FRI pop- 
ulation alone. As such understanding the little studied FRI sources 
and their evolution could be critical to deciphering this mechanism. 

The first significant attempt at determining the FRI high- 
redshift space density was carried out by Snellen & Best (2001) 
using the Hubble Deep Field and Flanking Fields (HDF+FF). Two 
z>l FRI galaxies were detected in this area, which calculations 
showed to be broadly consistent with an FRI space density en- 
hancement comparable to that of less luminous FRII galaxies at 
that redshift, and inconsistent (a probability of <1%) with a non- 
evolving FRI population. However, with only two detected FRIs 
the uncertainties in this result are clearly large. 

The area of sky used in the analysis of Snellen & Best (2001) 
was only large enough to give a first estimate of the high redshift 
space density of FRIs. This work, therefore, uses a deep, wide- 
field, Very Large Array (VLA) A-array survey an order of magni- 
tude larger than the HDF+HFF which will enable the space den- 
sity to be directly measured for the first time. Here we present the 
initial radio, optical and infra-red imaging, along with the spectro- 
scopic observations which form the basis of the work. The layout 
of the paper is as follows: in Sections[2]and[3]the sample is defined 
and the new radio observations taken are described; in Sections [4] 
and[5] the optical and infra-red imaging are presented and the host 
galaxies identified; Sections|6]and|7]describe the spectroscopic ob- 
servations of a subset of the sample. Finally, Section[8]outlines the 
redshift estimation methods used for the remainder of the sample 
and the conclusions of the paper can be found in Section[9] 



2 THE RADIO SAMPLE 

The survey was split over two fields - one in the constellation of 
Lynx at right ascension, a — 8/i45, declination, 5 — +44.6° 
(J2000) and one in Hercules at a = 17/i20, 8 = +49.9° (J2000). 
These fields were chosen because of the existence of previous low 
resolution radio observations by Windhorst et al. (1984), Oort & 
Windhorst (1985) and Oort & van Langevelde (1987). Addition- 
ally, the Hercules field has some previous optical and spectroscopic 
observations by Waddington et al. (2000). Alongside this, the Lynx 
field is also covered by the Sloan Digital Sky Survey (SDSS; York 
et al. 2000; Stoughton et al. 2002). 

The two fields were originally observed as part of the Leiden- 
Berkeley Deep Survey (LBDS) in which they were referred to as 
Lynx.2 and Hercules. 1. For the remainder of this paper they will be 
referred to as Lynx and Hercules respectively. It should be noted 
that the Lynx.2 field is unusual in that it does not contain any radio 
sources brighter than 6 mjy over nearly a square degree. However, 
the resulting underrepresentation in the radio source counts above 
this level should not affect the conclusions of this paper, as it is the 
faint end of the RLF that is being investigated here. This section 
outlines the previous work in the two fields. 



2.1 Sample definition and previous radio work 

The LBDS survey was constructed to provide photometry for faint 
galaxies and quasars, via multicolour plates obtained with the 4m 
Mayall Telescope at Kitt Peak (Kron, 1980; Koo & Kron 1982). 
Radio follow-up of nine of the LBDS fields (including Hercules 
and Lynx that are used here) was performed subsequently using the 
3km Westerbork Synthesis Radio Telescope (WSRT) at 1.4 GHz, 
at a resolution of 12.5" (Windhorst et al. 1984), reaching a rms 
noise level, at the field centre, of 0. 12-0.28 mjy. Their radio sample 
consists of 306 sources which satisfied the sample selection criteria 
of peak signal to noise (Sp /N) ^ 5a out to an attenuation factor, 
A(r), ^ 5 (for WSRT this corresponds to a radius of ^ 0.464° or 
«S 28'). 

The Hercules and Lynx fields were reobserved by Oort & 
van Langevelde (1987) and Oort & Windhorst (1985) respectively, 
again using the 3 km WSRT at 1.4 GHz with a 12.5" beam. These 
two sets of observations were a factor 2-3 deeper than the original 
Windhorst et al. (1984) ones, reaching a 5a flux limit of 0.45 mjy 
for Hercules and 0.30 mJy for Lynx at the pointing centre. 

The sample used in this work is a subset of the combined Her- 
cules and Lynx sources, as it is limited by the field of view size of 
the optical imaging described in Sj4j this is illustrated in FigureQ] A 
flux limit of 0.5 mjy was also imposed to remove the faintest, most 
poorly detected sources and provide a more uniform limiting flux 
density across the two fields. TableQ]gives the flux densities for the 
Hercules field (Oort & van Langevelde 1987) and Tabled gives the 
same for the Lynx field (Oort & Windhorst 1985) for the sources 
included in this work. Table [3] gives the 1.4 GHz flux densities of 
sources which were not covered by the optical observations, but 
were included in the subsequent infra-red imaging; consequently 
they are not part of the complete sample, but are included here 
for completeness. These, and other source parameters, were mea- 
sured by Oort and collaborators using an elliptical Gaussian fitting 
method in two dimensions. 



2.2 Sample completeness 

The selection criteria of the LBDS meant that source weighting 
was necessary to account for incompleteness in the starting sam- 
ple. This incompleteness arises from two factors: the attenuation of 
the WSRT primary beam and the resolution bias. The effect of the 
first of these factors, the decreasing sensitivity at increasing radial 
distance from the pointing centre, is to make the probability of de- 
tecting a source depend on where it is located in the map; a source 
which just satisfies the selection criteria at the centre would have 
been missed if it was located at the map edge. To correct for this 
each source was assigned a weight that was inversely proportional 
to the area over which it would have met the selection criteria and 
hence, have been included in the sample (Windhorst et al. 1984). 
The total sample area for the current work is limited by the size 
of the optical imaging, so the attenuation weights for the sources 
considered here were recalculated to account for this. 

The corrections for the incompleteness due to the resolution 
bias (the fact that a resolved source will be more difficult to de- 
tect than a point source of the same total flux) were found, by 
Windhorst et al. (1984), through detailed modelling of the source 
detection algorithms. It was subsequently found that sources with 
higher flux densities were also those with a large angular size (e.g. 
Windhorst et al. 1993) which, as Windhorst, Mathis & Neuschaefer 
(1990) showed, meant that the original weights were overestimated. 
Waddington et al. (2000) derived a new expression for the resolu- 



A sample ofmJy radio sources at 1.4 GHz in the Lynx and Hercules fields 3 




Figure 1. The distribution of the radio sources in the Lynx and Hercules fields; the labels correspond to the last 3 characters of the 
source names. The large circle in both plots represents the 0.464° radius of the previous WSRT observations of e.g. Windhorst et al. 
(1984); this is also approximately the 30' primary beam of the VLA at 1.4 GHz. The dotted line represents the area of the sample which 
is defined by the field of view of the optical imaging. Dots indicate the positions of the additional sources which were not included in 
the sample because they either fell below the flux density limit of 0.5 mJy, or they are not covered by the sample area. Open diamonds 
indicate sources not included in the complete sample, but which were included in the infra-red observations. 



tion bias that took this into account; it is that method which is used 
in the current calculations. 

These two re-derived correction factors, the attenuation 
weight and resolution weight, were then multiplied together to give 
the final weights. For further details of the source weighting meth- 
ods used in the LBDS see Windhorst et al. (1984) and Waddington 
et al. (2000). See Tables[]]and|2]for the weight calculated for each 
source. 



3 NEW RADIO OBSERVATIONS 
3.1 The observations 

The VLA data were taken with the array in A configuration on 22nd 
April 2002 for the Hercules field, and 15th February 2002 for the 
Lynx field. Both sets of observations were taken at a frequency of 
1.4 GHz (L-band), in four IF spectral line mode to enable wide field 
imaging, using 16 channels of width 781.25 kHz. The IFs were 
centred on 1391.3 and 1471.1 MHz observing dual polarization. 
Both the Hercules and the Lynx fields were observed for 8.5 hours 
each. 

The data were calibrated using the NRAO AIPS package. Be- 
cause of the non-coplanar geometry of the telescope array, in order 
to image the entire area of the primary beam, a three-dimensional 
Fourier transform is required. As this would be computationally 
impractical, the standard pseudo-three-dimensional Fourier trans- 
form technique, as incorporated in the task IMAGR, (Perley 1999) 
was adopted: this technique divides the field of view into numer- 
ous smaller facets, within each of which a two-dimensional Fourier 
transform provides a sufficiently good approximation. For each of 
the Hercules and Lynx fields, 256 by 256 pixel facets (with 0.35" 
per pixel) were centred on all the sources that were already known 
from the imaging of Windhorst, Oort and collaborators, discussed 
above. These sub-fields were cleaned and self-calibrated using the 
AIPS tasks IMAGR and CALIB. The self calibration consisted of 
multiple phase-only cycles followed by one final amplitude and 
phase calibration. The resulting maps have a resolution of 1.6" and 
reach a noise limit of 15 fiiy. 



3.2 Source detection and flux density measurement 

The Hercules sources of Oort & van Langevelde (1987) and Lynx 
sources of Oort & Windhorst (1985) were all detected in our VLA 
observations. This gives a sample of 81 sources, evenly spread over 
the two fields; the distribution on the sky can be seen in FigureQ] 
Flux densities were measured for these sources using the AIPS task 
imfit to fit a Gaussian, if they were pointlike, or the task tvstat to 
sum within a defined area, if they showed significant extension; 
the method used for each individual source is indicated in Tables 
Q]and[2] The values were then corrected for the attenuation of the 
VLA primary beam. For some objects, where previously one source 
was detected, these higher resolution observations have resolved 
it into more than one component each associated with a different 
host galaxy. In these cases, the sub-components are labelled a, b, 
etc. and the low resolution flux density has been assigned to each 
new component according to the A-array flux density ratio. These 
components are only retained in the sample if they remain above 
the 0.5 mJy flux density limit. 

Figure [2] shows a comparison between the WSRT flux den- 
sities measured by Oort et al. (1985; Lynx field) or Oort & van 
Langevelde (1987; Hercules field) and those measured for the 
new VLA A-array observations. The values are generally in good 
agreement for the compact sources; those that are different have 
•Soort > Si.4GHz which suggests that flux has been lost at the 
higher VLA resolution, and possibly indicates a resolved out FRI- 
type structure. The extended sources, however, are mainly underes- 
timated by Oort et al. (Soort < Si.4GHz) as a result of their ellipti- 
cal Gaussian measuring method which misses any extended flux. It 
should be noted that a subset of the two fields (36% of the sample) 
were also observed with the VLA by Oort et al. (1987). These mea- 
surements have not been used for flux density comparisons though, 
due to the small number of sources included. 

The positions of the detected sources are given in Tables[T]and 
|2]along with the measured flux densities and primary beam correc- 
tion factors, Cpb, (i.e. Scor = Cpb Sivieas) used. The correspond- 
ing radio contour images can be found in Appendix iBl 
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Table 1. The Hercules radio source positions from the VLA A-array observations along with the A-array and Oort 
et al. 1987 1.4 GHz primary beam corrected flux densities, source weights and primary beam correction factor, Cpb 
(see text for full details). An T in the final column indicates an imfit measured flux density; a T indicates a tvstat 
measurement. 



Hercules 


Name 


RA/DEC (J2000) 


Srtnrt fmJv) 


S^ inn, (mjy) 


Weight 


CpB 


Measure 


53w052 


17 18 34.14 


a r\ c n C"i n 

49 58 53.0 


8.00 ± 0.34 


9.31 ± 1.17 


1.00 


2.26 


I 


53w054a 


1 "7 10 Al i n 

17 lo 47.30 


a r\ ac a n A 

4y 45 4y.u 


2.07 ±0.19 


2.14 ±0.35 


1.00 


2.34 


I 


53w054b 


1 "7 1 o a n CV~I 

17 lo 4v.y7 


An A £. 111 

49 4o 12.2 


2.08 ±0.19 


2.44 ± 0.32 


1.00 


2.20 


I 


53w057 


17 19 07.29 


AC\ AC A A O 

49 45 44.8 


1.96 ±0.14 


1.95 ±0.21 


1.00 


1.82 


I 


53w059 


\ i 1 n in 10 

1 7 19 20. lo 


c n c\f\ 1 1 i 

50 00 21.2 


19.40 ± 1.0 


23.81 ± 1.21 


1.00 


1.34 


T 


53w061 


17 19 27.34 


49 44 01.9 


4.76 ± 0.43 


1.44 ±0.18 


1.02 


1.69 


I 


53w062 


17 19 31.9 j 


Af~\ Cd 1 

49 59 06.2 


0.73 ±0.10 


1.08 ±0.07 


1.42 


1.19 


I 


53w065 


1 'H in A f\ l\C 

17 19 40.05 


49 57 39.2 


5.54 ± 0.20 


5.89 ±0.14 


1.00 


1.11 


I 


53w066 


17 19 42.96 


50 01 03.9 


4.27 ±0.17 


4.53 ±0.15 


1.00 


1.18 


I 


53w067 


17 19 51.27 


50 10 58.7 


21.9 ±0.90 


36.68 ± 3.97 


1.00 


2.15 


T 


53w069 


17 20 02.52 


49 44 51.0 


3.82 ±0.17 


5.25 ±0.31 


1.00 


1.36 


T 


53w070 


17 20 06.07 


50 06 01.7 


2.56 ±0.14 


2.61 ±0.17 


1.00 


1.39 


I 


53w075 


1 t in a*~\ n 

17 20 42.37 


AC\ A^i AC\ 1 

49 43 49.1 


96.8 ± 3.3 


99.82 ± 6.83 


1.00 


1.51 


I 


53w076 


17 20 55.82 


49 41 02.2 


1.94 ±0.17 


6.93 ± 0.92 


1.00 


2.21 


T 


53w077 


17 21 01.32 


Af~\ A O 1 A f\ 

49 48 34.0 


6.51 ±0.39 


18.11 ± 1.01 


1.00 


1.31 


T 


53w078 


17 21 18.17 


50 03 35.2 


0.74 ±0.12 


1.84 ±0.20 


1.40 


1.62 


T 


53w079 


In 11 *"7C 

17 21 22.75 


rn in ti f\ 

50 10 31.0 


11.7 ± 0.5 


11.1 ± 1.68 


1.00 


2.87 


I 


53w080 


17 21 37.48 


Af\ C C 1 £. O 

49 55 36.8 


25.9 ± 0.9 


31.11 ± 2.42 


1.00 


1.63 


T 


53w081 


17 21 37.86 


AC\ C ~T C~l f 

49 57 57.6 


12.1 ± 0.5 


12.93 ± 1.08 


1.00 


1.68 


I 


53w082 


17 21 37.64 


50 08 27.4 


2.50 ±0.19 


2.97 ± 0.47 


1.00 


2.86 


I 


53w083 


17 21 48.95 


50 02 39.7 


5.01 ± 0.25 


5.06 ± 0.64 


1.00 


2.28 


I 


53w084 


17 21 50.43 


49 48 30.5 


0.68 ±0.12 


1.01 ±0.19 


1.51 


2.53 


I 


53w085 


17 21 52.48 


49 54 34.1 


4.52 ± 0.22 


4.94 ± 0.66 


1.00 


2.02 


I 


53w086a 


17 21 56.42 


49 53 39.8 


1.62 ±0.30 


4.06 ± 0.54 


1.09 


2.17 


T 


53w086b 


17 21 57.65 


49 53 33.8 


2.44 ± 0.30 


6.13 ±0.74 


1.00 


2.22 


T 


53w087 


17 21 59.10 


50 08 42.9 


5.58 ±0.35 


14.35 ± 2.23 


1.00 


4.23 


T 


53w088 


17 21 58.90 


50 11 52.7 


14.1 ± 0.7 


14.52 ± 2.92 


1.00 


6.07 


I 


53w089 


17 22 01.05 


50 06 54.7 


3.04 ± 0.26 


3.58 ± 0.62 


1.00 


3.71 


T 


66w009a 


17 18 32.76 


49 55 53.4 


1.14 ±0.21 


1.50 ±0.22 


1.23 


2.22 


I 


66w009b 


17 18 33.73 


49 56 03.2 


0.70 ± 0.21 


0.91 ±0.16 


5.69 


2.19 


I 


66w014 


17 18 53.51 


49 52 39.1 


3.34 ±0.51 


0.60 ± 0.09 


1.18 


1.66 


I 


66w027 


17 19 52.11 


50 02 12.7 


0.57 ±0.11 


0.67 ±0.13 


3.75 


1.19 


I 


66w031 


17 20 06.87 


49 43 57.0 


0.76 ±0.14 


0.97 ±0.12 


2.32 


1.43 


I 


66w035 


17 20 12.32 


49 57 09.7 


0.63 ± 0.09 


0.71 ±0.06 


1.77 


1.01 


I 


66w036 


17 20 21.46 


49 46 58.3 


0.78 ±0.11 


3.70 ± 0.29 


1.57 


1.20 


T 


66w042 


17 20 52.59 


49 42 52.4 


0.78 ±0.14 


1.99 ±0.26 


1.61 


1.70 


T 


66w047 


17 21 05.43 


49 56 56.0 


0.60 ±0.10 


1.16 ± 0.10 


2.32 


1.20 


T 


66w049 


17 21 11.25 


49 58 32.4 


1.38 ±0.27 


2.17 ±0.22 


3.09 


1.28 


I 


66w058 


17 21 48.23 


49 47 07.3 


1.89 ±0.16 


1.72 ±0.24 


1.01 


2.33 


I 



4 OPTICAL AND INFRA-RED IMAGING 



The FRI space density calculation depends on determining the red- 
shifts of the radio galaxies in the survey. To obtain all of these spec- 
troscopically would have been very time consuming, so the aim 
was to combine photometric redshift estimates with spectroscopic 
follow-up of the best high redshift FRI candidates. Optical obser- 
vations were carried out using the Wide Field Camera (WFC) on 
the 2.5 m Isaac Newton Telescope (INT) in La Palma, and a sub- 
sample of sources were also observed using the UKIRT Fast Track 
Imager (UFTI) on UKIRT, the 3.8 m UK Infra-red Telescope lo- 
cated in Hawaii. These two sets of observations are described in 
this section. 



4.1 INT observations and data reduction 

The WFC consists of 4 thinned EEV 2kx4k CCDs with a pixel size 
of 13.5 /im, resulting in a scale of 0.33" /pixel and a combined field 
of view of ~ 34 x 34 arcmin 2 . This large field of view makes the 
WFC an ideal instrument for observing the two fields which are of 
comparable size. 

The main WFC observations were split over two separate 
runs in April 2003 and April 2004. Unfortunately these were both 
largely weathered out, so observations through two filters only were 
obtained - sloan r 1 and (', Several exposures of 300s or 600s were 
taken, and the telescope was offset by 30" after every third obser- 
vation to fill the gaps between the CCDs and cover the whole field. 
Full details of the observations can be found in Table|4] 

The April 2003 run took data on one night only; the 7th. Con- 
ditions were non-photometric, so these data have only been used to 
determine optical identifications for the radio objects. For the April 
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Table 2. The Lynx radio source positions from the VLA A-array observations along with the A-array and Oort et 
al. 1985 1.4 GHz primary beam corrected flux densities, source weights and primary beam correction factor, Cpb 
(see text for full details). An T in the final column indicates an imfit measured flux density; a 'T' indicates a tvstat 
measurement. 



Lynx 


Name 


RA/DEC (J2000) 


Sn^rt fmJv) 

"Uon v J j t 


S^ An u-7 (mJv) 


Weight 


Cp B 


Measure 


55wll6 


no An a f\ to 

08 43 40.72 


a A on o A i 

44 39 24.7 


1.36 ±0.12 


1.52 ±0.25 


1.00 


2.22 


I 


55wll8 


no An A f O/T 

08 43 46.86 


A a o c 1 r\ i 

44 35 49.7 


0.91 ± 0.09 


0.74 ±0.11 


1.03 


1.92 


I 


55wl20 


AO Al on 

08 43 52.89 


A A O A Of\ f\ 

44 24 29.0 


1.83 ±0.16 


1.67 ±0.29 


1.00 


2.68 


I 


55wl21 


no a A n a f\/~ 

08 44 04.06 


A A o 1 1 n a 

44 31 19.4 


1.21 ±0.09 


1.04 ±0.11 


1.00 


1.60 


I 


55wl22 


no A A 1 O O O 

08 44 12.33 


A A ^ 1 1 a n 

44 31 14. 9 


0.56 ± 0.08 


0.66 ±0.12 


1.27 


1.45 


I 


55wl23 


no A a 1 A c a 

08 44 14.54 


a a o c nn o 

44 35 00.2 


2.01 ±0.10 


1.17 ±0.08 


1.00 


1.33 


I 


55wl24 


no A a a A no 

08 44 14.93 


A A O O CO O 

44 38 52.2 


4.67 ±0.17 


2.79 ±0.16 


1.00 


1.35 


I 


55wl27 


no A A ot c c 

08 44 27.55 


A A An (VI A 

44 43 07.4 


1.81 ±0.10 


1.64 ±0.11 


1.00 


1.36 


I 


55wl28 


no A A O O Af 

08 44 33.05 


A A C f\ 1 C O 

44 50 15. 3 


3.34 ±0.18 


4.77 ± 0.54 


1.00 


2.05 


T 


55wl31 


no A A n r r i 

08 44 35.51 


A A A f f\ A 1 

44 46 04. 1 


1.01 ±0.10 


0.74 ±0.11 


1.01 


1.48 


T 


55wl32 


no A A O T -i o 

08 44 37.12 


A A C f\ O A 1 

44 50 34.7 


1.10 ±0.11 


1.66 ±0.23 


1.01 


2.05 


T 


55wl33 


nO A A ni n A 

08 44 37.24 


A A O S~ AA A 

44 26 00.4 


2.20 ±0.11 


2.25 ±0.16 


1.00 


1.47 


I 


55wl35 


no A A A 1 1 n 

08 44 41.10 


A A O 1 O "7 O" 

44 21 37.7 


2.60 ±0.14 


3.86 ± 0.44 


1.00 


1.98 


T 


55wl36 


08 44 45.14 


44 32 23.9 


1.02 ±0.07 


0.92 ± 0.08 


1.00 


1.10 


T 


55wl37 


no A a a f nn 

08 44 46.90 


A A A A n^l n 

44 44 37.9 


1.60 ±0.09 


1.66 ±0.11 


1.00 


1.29 


T 


55wl38 


08 44 54.51 


44 46 22.0 


1.82 ±0.10 


1.99 ±0.15 


1.00 


1.37 


I 


55wl40 


no ac f\£. f\f 

08 45 06.06 


A A Af\ A 1 O 

44 40 41. 2 


0.79 ± 0.08 


0.55 ± 0.06 


1.25 


1.06 


I 


55wl41 


no ac no on 

08 45 03.29 


A A OO 1 C 1 

44 28 15.1 


0.87 ± 0.07 


0.43 ± 0.06 


1.01 


1.19 


I 


55wl43a 


no ac n c An 

08 45 05.49 


A a o c ,1 c n 

44 25 45.0 


2.41 ± 0.11 


2.19 ±0.13 


1.00 


1.34 


I 


55wl43b 


no AC f\ A or 

08 45 04.25 


A A O C CO O 

44 25 53.3 


0.57 ± 0.09 


0.33 ± 0.06 


1.58 


1.33 


I 


55wl47 


no Ac oo oo 

08 45 23.83 


I A C f\ n A s 

44 50 24.6 


1.72 ±0.11 


1.97 ±0.19 


1.00 


1.82 


I 


55wl49 


08 45 27.17 


44 55 25.9 


7.10 ± 0.32 


7.82 ± 1.11 


1.00 


3.20 




55wl50 


no /icon a "i 

08 45 29.47 


A a cc\ n*i A 

44 50 37.4 


0.95 ± 0.10 


0.63 ± 0.10 


1.02 


1.88 


I 


55wl54 


08 45 41.30 


44 40 11.9 


12.1 ± 0.40 


13.71 ± 0.40 


1.00 


1.13 




55wl55 


08 45 46.89 


44 25 11.6 


1.83 ± 0.10 


1.70 ± 0.14 


1.00 


1.55 


I 


55wl56 


08 45 50.92 


44 39 51.5 


4.14 ±0.16 


4.78 ± 0.21 


1.00 


1.19 




55wl57 


08 46 04.44 


44 45 52.7 


1.37 ±0.10 


1.24 ± 0.12 


1.00 


1.68 




55wl59a 


08 46 06.67 


44 51 27.5 


6.70 ± 0.29 


6.49 ± 0.82 


1.00 


2.69 




55wl59b 


08 46 06.82 


44 50 54.1 


0.75 ± 0.13 


1.00 ±0.19 


1.00 


2.55 




55wl60 


08 46 08.50 


44 36 47.1 


0.94 ± 0.08 


0.81 ± 0.07 


1.01 


1.32 




55wl61 


08 46 27.32 


44 29 56.9 


1.34 ±0.14 


1.25 ± 0.15 


1.02 


1.87 




55wl65a 


08 46 34.76 


44 41 39.2 


18.12± 0.54 


18.88 ± 1.54 


1.00 


2.06 




55wl65b 


08 46 33.37 


44 41 24.4 


0.78 ± 0.40 


0.92 ± 0.14 


1.47 


1.99 




55wl66 


08 46 36.02 


44 30 53.5 


2.46 ±0.14 


2.31 ± 0.26 


1.00 


2.07 




60w016 


08 44 03.58 


44 38 10.2 


0.62 ± 0.08 


0.88 ±0.13 


1.14 


1.52 




60w024 


08 44 17.83 


44 35 36.9 


0.51 ±0.09 


0.37 ± 0.05 


2.06 


1.29 




60w032 


08 44 33.69 


44 46 13.0 


0.54 ± 0.09 


0.46 ± 0.08 


1.58 


1.51 




60w039 


08 44 42.50 


44 45 32.5 


0.65 ± 0.09 


0.72 ±0.16 


1.26 


1.38 




60w055 


08 45 14.00 


44 53 08.7 


0.51 ± 0.08 


0.62 ±0.13 


1.42 


2.34 




60w067 


08 45 40.47 


44 23 20.1 


0.56 ± 0.09 


0.69 ±0.15 


1.31 


1.70 


T 


60w071 


08 46 00.34 


44 43 22.1 


0.50 ± 0.08 


0.60 ± 0.07 


1.44 


1.42 


I 


60w084 


08 46 39.86 


44 33 44.5 


0.85 ± 0.17 


1.80 ±0.37 


1.38 


2.08 


T 



Table 3. The Hercules and Lynx field sources which did not fall within the optical field but were included in the 
infra-red observations. They are given a weight of 0.00 as they are not part of the complete sample. They are included 
here for completeness. For full details see text. An T in the final column indicates an imfit measured flux density; a 
T indicates a tvstat measurement. 



Name 


RA/DEC (J2000) 


^oort (mJy) 


Si.4GHz (mJy) 


Weight 


Cpb 


Measure 


53w091 


17 22 32.73 


50 06 01.9 


22.6 ± 1.1 


37.93 ± 6.62 


0.00 


7.43 


T 


55wll9 


08 43 47.98 


44 50 41.4 


1.78 ± 0.16 


1.78 ± 0.30 


0.00 


3.86 


I 


55wl25 


08 44 15.28 


44 11 16.2 


22.2 ± 1.20 


12.15 ±2.30 


0.00 


12.01 


I 


55wl26 


08 44 20.56 


44 58 05.0 


3.20 ± 0.24 


4.19 ±0.76 


0.00 


6.37 


I 
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Table 4. Summary of INT observations 



Field Band Observation Date Exposure Time Photometric? Seeing (") 



Hercules r 07/04/03 24x300s No 2.5 

i " 10x300s No 1.9 

r 15/04/04 15x600s Yes 1.5 

i " 15x300s Yes 1.5 

Lynx i 07/04/03 9x300s No 2.1 

r 15/04/04 6x300s Yes 1.5 

i 06/01/05 1x300s Yes 3.0 



1 10 100 

So* / r"Jy 



Figure 2. A comparison of the VLA A-array (Si.4GHz) ar| d Oort et al. 
(1985) or Oort & van Langevelde (1987) (Soort) flux densities for both 
fields. Extended, tvstat measured sources are shown as triangles and com- 
pact, imfit measured sources are shown as diamonds. 

2004 run useful observations were only taken on the 15th. Standard 
star fields were observed throughout the photometric night. 

One further observation of the Lynx field of 300s, in was 
obtained on 6th January 2005. This was needed because of the lack 
of photometric i-band data in this field. Although the night was 
photometric the seeing was very poor (~ 3"). Standard star fields 
were again observed throughout the night. This image was only 
used to photometrically calibrate the Lynx field; the identifications 
were done with the previous r and i-band images. 

All the images were processed using the IRAF software pack- 
age. Bias frames, taken at the beginning of each night, were aver- 
aged together to make a master-bias for each detector which was 
then subtracted from the remaining data. For the majority of the ob- 
servations, a flat-field was made for the four detectors, in each field 
and filter, by median combining the separate science frames and re- 
jecting pixels according to the readnoise and gain of the CCD. The 
exception to this was the 15th April observations where twilight 
flats were taken. Next the individual science frames were divided 
by the corresponding sky-flat, which had been normalised using its 
pixel mean. The frames were registered using ~10 stars and shifted 
and the final images were then median combined and clipped us- 
ing the CCD noise properties as before. The offsets between frames 
were sufficiently small that it was not necessary to account for dis- 
tortions in this process. 

The astrometric calibration for the INT data was complicated 
by the distortion of the WFC across the 4 CCDs (Taylor, 2000). To 
correct for this the USNcQBl.O catalogue (Monet et al. 2003) was 

1 United States Naval Observatory 



used to create a catalogue of reference stars, which was then used 
by the starlink program astrom to calculate a distortion-corrected 
calibration. For the April 2004 observations this astrometric cali- 
bration was only applied to the r-band data, as it was felt that no 
further radio host-galaxy identifications would be obtained from 
considering the significantly shallower /-band data also. However, 
the i-band images were tweaked locally to each r-band detection to 
ensure that the images lined up. The calculated errors were ~ 0.3" 
for all CCDs in both fields. 



4.2 UKIRT observations and data reduction 

In contrast, UFTI consists of one 1024x1024 HgCdTe array with a 
plate scale of 0.09l"/pixel. This results in a field of view of 92" 
which is significantly smaller than that of the WFC. This meant 
that it could only be used to obtain images of, mainly, individual 
sources rather than the complete-field observations done with the 
INT. 

The UKIRT UFTI observations were done in a combination 
of service and visitor mode, spread over the period July 2004 to 
January 2005. All the observations were done using the AT-band 
filter; details of the observations, along with the sources observed 
can be found in Table|5] The sources selected for these observations 
were those with a faint optical detection or no optical detection at 
all. 

The sources were observed using a 9-point dither pattern, with 
offsets of 10", and an exposure time of 60s per dither position. In 
general this was repeated 4 times resulting in a total of 36 expo- 
sures. The exceptions to this were 53w054, where the observation 
had to be re-started due to high humidity and 55wl25, where soft- 
ware problems meant the observation had to be halted after 27 ex- 
posures (3 repeats of the dither pattern). Also, because some obser- 
vations were done in service, two sources which were not detected 
after 36 exposures were able to be re-observed for a further 18 ex- 
posures at a later date. The observations of 55wl28 and 55wl32 
were also repeated since the originals were taken at a very high air- 
mass which resulted in significant elongation of the objects in the 
field. 

Appropriate standard stars (FS 125 for Lynx and FS27 for Her- 
cules) were observed multiple times throughout the night if multi- 
ple targets were also observed, but only observed once on other 
nights. All nights were photometric. 

The infra-red data reduction method is similar to that already 
described for the optical data; again the IRAF software was used 
to process the images, but each individual source observation was 
reduced independently. The first step was the subtraction of the ap- 
propriate dark frame from each image. Flat-fields were then made 
for each source by median combining the first 9 observations only, 
rejecting pixels according to the noise properties of the detector, to 
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Table 5. Summary of UKIRT observations. 



Table 6. The calibration co-efficients for the two observations 



Field 



Lynx 



Date 


Target 
Source 


Exposure 
Time 


Seeing 


25/07/04 


53w054 


54x60s 


0.7 


28/07/04 


53w084 


36x60s 


0.5 




53w087 


36x60s 




22/08/04 


53w089 


36x60s 


0.9 




66w031 


36x60s 




12/09/04 


66w009 


36x60s 


0.8 


1 1/09/04 


53w091 


36x60s 


0.7 


14/09/04 


66w035 


36x60s 


1.1 


20/09/04 


66w036 


36x60s 


0.9 


15/01/05 


55wll9 


36x60s 


0.9 


21/01/05 


55wl28 


36x60s 


0.7 




55wl32 


36x60s 






55wl20 


36x60s 






55wl25 


27x60s 






55wl26 


36x60s 






55wl33 


36x60s 




23/01/05 


55wl35 


36x60s 


0.7 




55wl38 


36x60s 






55wl47 


36x60s 




24/01/05 


55wl55 


36x60s 


0.5 




55wl36 


36x60s 






55wl28 


36x60s 






55wl32 


36x60s 




16/02/05 


55wl33 


18x60s 


0.9 


17/02/05 


55wl21 


36x60s 


1.4 




55wl23 


36x60s 






55wl56 


36x60s 






55wl43 


36x60s 





minimise the effects of sky variability over the full length of the ex- 
posure. The images were then divided by the flat-field which had 
been normalised using its median pixel value. 

The final steps in the reduction process - sky-subtraction, cos- 
mic ray removal and image combining - were done using the IRAF 
package dimsum, by P. Eisenhardt, M. Dickinson and S.A. Stan- 
ford, and, in particular, the task reduce. The images were again 
registered using, on average, 10 stars. For the source with 2 sep- 
arate observations (55wl33) the sky-subtraction was done for the 
two nights separately, but all the images were then registered and 
combined together to produce a single final image. 

The astrometric calibrations for the images were done where 
possible using the INT images as references. In cases with no INT 
overlap, rough astrometry was derived from the telescope pointing 
position and then the RA and DEC positions were improved using 
the one or two stars from the USNO catalogue available in the small 
fields. 



4.3 Aperture photometry and source identification 

This section outlines the steps taken to identify the host galaxies of 
the sample sources and the subsequent magnitude measurements of 
these objects. 



4.3. 1 Optical standard star calibration 

The standard star observations were reduced in the same way as the 
science observations. The fields used were from the Landolt Faint 
Equatorial Standards catalogue (Landolt, 1992), each containing 



Date Filter 


n (mag/airmass) 


m zp t 


April 2004 i 


-0.03 


24.24 ± 0.05 


r 


-0.07 


24.68 ± 0.05 


January 2005 i 


-0.01 


24.31 ±0.05 



an average of 10 standard stars. The April 2004 standard field was 
SA104 in r-band only and SA107 in r and i, whereas the January 
2005 standard fields were SA104 and SA98. The April 2003 obser- 
vations were not photometric. Object counts were measured using 
a 5" radius aperture for all standards using the gaia package; the 
only exception to this was for the January 2005 standards where 
high seeing meant a larger aperture (15" radius) was needed, and 
stars with near neighbours were ignored to minimise errors. 

The Landolt Faint Equatorial Standards are only available in 
the Johnson-Rron-Cousins photometric system, and are given in 
Vega magnitudes. Therefore the apparent magnitudes for the Lan- 
dolt standard stars needed to be transformed to the Sloan photomet- 
ric system. This transformation was done using the following two 
relations from Smith et al., (2002): 



V - 0.84(1/ - R) +0.13 
-% = 1.00CR-7) -0.21 



(1) 
(2) 



In the above equations lowercase letters indicate the Sloan photo- 
metric system and uppercase letters the Johnson system. It should 
be noted that this also transforms the magnitudes to the AB- 
magnitude system. 

Once the transformations had been applied, the calibration co- 
efficients, zeropoint magnitude, m zpt , and extinction co-efficient, 
k, were determined for the photometry. These are summarised in 
Tableland are in good agreement with previously published val- 
ues for INT extinction. 



4.3.2 Infra-red standard star calibration 

The UKIRT standard star observations were reduced in a similar 
way to the science images, the one difference being that all 5 ob- 
servations for each standard were used to make the flat-field image. 
In contrast to the INT Landolt standard star fields, which contained 
many stars, only one star was used for Hercules (FS27) and one for 
Lynx (FS125). These were both from the UKIRT Faint JHK Stan- 
dards catalogue and are given in Vega magnitudes (Casali, 1992). 
The aperture size used to measure the standards was 2.5" radius. 

As there is only one star per field and, in many cases, just one 
standard star observation per night, the zeropoint magnitude and 
extinction co-efficient could not both be determined. The extinc- 
tion co-efficient, n, was therefore taken to be 0.05 mags/airmass, 
the published value for UFTI (Leggett, 2005). The zeropoint mag- 
nitudes for each observing night were then calculated using the 
published m app for each standard; the values found all lay in the 
range 22.35-22.40 and are in good agreement with previous values 
given for UKIRT. For nights where more than one source and hence 
more than one standard, were observed the mean value for m zp t 
was used; the uncertainty of ±0.02 on this zeropoint was incorpo- 
rated into the error estimates of the source magnitudes as described 
in 34331 
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4.3.3 Identifications and magnitudes 

The source host-galaxy positions were found by overlaying the 
VLA A-array radio contour maps with the optical and infra-red 
data. Figures B 1 and B2 show the radio/optical and, where appro- 
priate, radio/infra-red overlays resulting from the UKIRT and INT 
April 2004 observations. Radio contour maps for sources with no 
host galaxy identification are shown in Figure 3. The corresponding 
host galaxy positions are given in Tables [8] and [9] 

The aperture photometry of all the sources was then done, 
again using the gaia package. The counts received for each source 
in r, i, and K-band (if available), were measured in 4 different sized 
apertures - 1.5", 2.5", 4.0" and 8.0" radius - and then, depend- 
ing on the extent of the source, the measurement from the most 
appropriate aperture was selected and used thereafter. The aperture 
chosen for each source was the same in the three bands to enable 
colours to be accurately determined. Sky-subtraction was achieved 
either using an annulus round the object or, in cases where this 
could not be done because of the proximity of other objects, using 
a sky-aperture placed nearby. 

The magnitudes were then calibrated using the appropriate 
values of n and m zpt determined for the optical and infra-red stan- 
dard stars; these can be found in Tables [8] and [9] and are in rea- 
sonable agreement with previously published results. The AT-band 
magnitudes for sources observed with UKIRT but not included in 
the complete sample are given in Table|7] 

4.3.4 Aperture corrections 

The next step was to correct all the calculated source apparent mag- 
nitudes to a metric aperture of 63.9 kpc diameter, thus allowing ac- 
curate comparisons to be made between sources at all redshifts. The 
63.9 kpc aperture, corresponding to an aperture of ~8" at z = 1 has 
become a standard metric size following previous work by Eales et 
al. (1997) and others. 

At low redshift (z < 0.6) this correction is carried out us- 
ing the curve of growth for elliptical galaxies tabulated by Sandage 
(1972); this method assumes that the hosts of radio galaxies are all 
giant ellipticals and that they share the same intensity profile. This 
assumption is a good approximation at low redshift, but is not valid 
for higher redshift radio galaxies which can exhibit very different 
structure. For radio galaxies located at z > 0.6 therefore the mea- 
sured emission within an aperture of radius r was assumed to be 
proportional to r a where a — 0.35 (Eales et al., 1997). 

These aperture correction methods obviously depend on the 
redshifts for the sources being known; only a small fraction of the 
sample satisfied this condition. For the remaining objects redshifts 
were estimated iteratively using the K-z and r-z relations, as de- 
scribed in Sj8] The calculated magnitude corrections for the r, i and 
A"-band magnitudes can be found in Tables [8] and [9] The aperture 
corrections range from +0.15 to -0.52 magnitudes and are typi- 
cally negative, with an average correction of -0.3 mag to account 
for missing flux. 

4.3.5 Magnitude error 

The magnitude error was determined by combining in quadrature 
four potential sources of error: (i) the error on the received counts as 
determined by Poisson error on the incoming photons, (ii) the error 
in the determination of m zpt , (iii) the error due to the subtraction of 
the sky background, found by taking the standard deviation of 10 
apertures placed randomly on empty regions in the fields and (iv) 



Table 7. The host galaxy positions and A'-band magnitudes for the sources 
not included in the complete sample. The radius (in ") of the aperture used 
for photometry is given in brackets and 3a limits are given for undetected 
sources. The corresponding radio positions can be found in TablefJ] 



Name 


RA 


DEC 


K 


K (63.9 kpc) 




(J2000) 


(J2000) 






53w091 


17 22 32.73 


50 06 01.9 


18.40 ±0.15 (2.5) 


18.25 ±0.17 


55wll9 


08 43 47.98 


44 50 41.4 


19.54 ± 0.37 (2.5) 


19.36 ±0.38 


55wl25 


08 44 15.25 


44 11 16.7 


17.44 ± 0.06 (2.5) 


17.26 ±0.11 


55wl26 






>19.85 





the error in the aperture correction which is taken as 50% of the 
correction value. If the source is bright (i) dominates; (iii) is most 
important for the infra-red observations where the background is 
very high. 



5 IMAGING RESULTS 

The April 2003 optical data resulted in an identification fraction of 
53% and 63% for the Lynx and Hercules fields respectively. These 
numbers rose to 76% and 87% with the inclusion of the optical 
data from April 2004. 80% of the Hercules sources observed in 
the infra-red were identified, compared with 57% of the sources 
observed in the Lynx field. Combining these figures gave an iden- 
tification fraction for the Lynx field of 83%, and 90% for Hercules. 

In total, out of the complete sample, 4 radio sources in the Her- 
cules field and 7 radio sources in the Lynx field remain unidentified 
after the r, i' and K-band observations. The observations reached 
optical 3<t limiting magnitudes of r' < 25.17 mag and i' < 23.76 
mag for Hercules, and r' < 24.38 mag and i' < 23.46 mag for 
Lynx; the infra-red 3cr limiting magnitudes were K < 19.85 mag, 
K < 19.98 mag and K < 20.16 mag for the 36x60s, 54x60s and 
72x60s observations respectively. Notes on individual sources can 
be found in Appendix lAll 

5.1 Colours and magnitude distribution 

The (r — i) and (r — K) colour-magnitude diagrams for both Lynx 
and Hercules sources are shown in Figure(3] The (r — i) plot sug- 
gests a slightly greater range in the colours of the sources in the 
Lynx field compared to the Hercules field, but this is likely to be 
the result of the relative shallowness of the Lynx i-band data, which 
would introduce a bias against the fainter bluer sources in this field, 
rather than a result of the aforementioned radio source count under- 
representation. The (r — i) colours for both fields though show that 
the majority of the sources are hosted by red galaxies as expected. 

The Lynx field sources in the (r — K) plot have a similar 
colour distribution to the Hercules sources, but the small number 
of sources in this diagram make a comparison difficult. The small 
numbers and poor population are the result of the selection crite- 
ria used for the K-band observations; sources with either an unde- 
tected or faint r-band detection were those chosen. 

Also shown on Figure are two likely models for an L* type 
galaxy (assuming M* = —22.0 (Blanton et al., 2001)), calculated 
using the GALAXEV code (Bruzual & Chariot 2003): passive evo- 
lution following an initial burst of star formation at z = 5, and 
exponentially declining star formation, with an e-folding time of 1 
Gyr, beginning at z = 10. The (r — i) vs i model lines track the 
data well out to !~20 (corresponding to z ~ 0.6), beyond which 
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Table 8. The host galaxy positions and magnitudes for the Hercules field. The radius (in ") of the aperture used for photometry is given in brackets. Sources 
which were unmeasurable due to the presence of a nearby bright object are labelled with a * and 3<r limits are given for undetected sources. The corresponding 
radio positions can be found in Table[T] 



Hercules 



Name 


T> A 
KA 

yj ZUUU ) 


yj ZL/UU ) 


r 


r \p j. v Kpcj 


i 


i (f^X O l-j-i/-"! 

1 [Os>.y KpCJ 


IV 


iv (^oj.y Kpcj 


53w052 


17 18 34.07 


49 58 50.2 


21.31 ± 0.05 (4) 


21.19 ±0.08 


20.86 ± 0.07 (4) 


20.74 ± 0.09 


_ 


_ 


53w054a 


17 18 47.30 


49 45 49.0 


23.74 ±0.14 (2.5) 


23.58 ±0.16 


23.62 ± 0.26 (2.5) 


23.46 ± 0.27 


18.32 ±0.13 (2.5) 


18.17 ±0.15 


53w054b 


17 18 49.97 


49 46 12.2 


>25.17 


_ 


>23.76 


_ 


19.95 ± 0.59 (2.5) 


19.75 ± 0.60 


53w057 


17 19 07.29 


49 45 44.8 


24.69 ±0.31 (2.5) 


24.53 ± 0.32 


>23.76 


_ 


_ 


_ 


53w059 


17 19 20.26 


50 00 19.6 


24.32 ± 0.22 (2.5) 


24.17 ±0.23 


>23.76 


_ 


_ 


_ 


53w061 


17 19 27.34 


49 43 59.7 


21.13 ± 0.05 (4) 


21.13 ± 0.05 


20.77 ± 0.07 (4) 


20.77 ± 0.07 


_ 


_ 


53w062 


17 19 32.07 


49 59 06.8 


21.91 ±0.06(2.5) 


21.67 ±0.14 


21.04 ±0.06 (2.5) 


20.80 ±0.14 


_ 


_ 


53w065 


17 19 40.07 


49 57 40.8 


23.00 ± 0.08 (2.5) 


22.84 ± 0.11 


23.31 ± 0.20 (2.5) 


23.14 ±0.22 


_ 


_ 


53w066 


_ 




>25.17 


_ 


>23.76 


_ 


_ 


_ 


53w067 


17 19 51.27 


50 10 58.5 


22.15 ±0.06 (2.5) 


21.94 ± 0.12 


21.43 ±0.06 (2.5) 


21.22 ±0.12 


_ 


_ 


53w069 


17 20 02.54 


49 44 51.0 


25.12 ±0.46 (2.5) 


24.97 ± 0.47 


>23.76 


_ 




_ 


53w070 


17 20 06.07 


50 06 01.7 


22.20 ± 0.06 (2.5) 


22.05 ± 0.10 


21.37 ±0.06 (2.5) 


21.21 ± 0.10 


_ 


_ 


53w075 


17 20 42.36 


49 43 49.2 


21.10 ±0.05 (4) 


21.12 ±0.05 


20.67 ± 0.06 (4) 


20.69 ± 0.06 


_ 


_ 


53w076 


17 20 55.78 


49 41 03.1 


19.57 ± 0.05 (4) 


19.41 ± 0.10 


18.91 ± 0.05 (4) 


18.75 ± 0.10 


_ 


_ 


53w077 


17 21 01.32 


49 48 34.1 


21.71 ± 0.05 (4) 


21.69 ±0.05 


20.82 ± 0.07 (4) 


20.80 ± 0.07 


_ 


_ 


53w078 


17 21 18.17 


50 03 34.9 


18.28 ± 0.05 (8) 


18.29 ± 0.05 


17.54 ± 0.05 (8) 


17.54 ± 0.05 


_ 


_ 


53w079 


17 21 22.62 


50 10 31.2 


20.62 ± 0.05 (4) 


20.54 ± 0.07 


19.71 ±0.05 (4) 


19.62 ± 0.07 


_ 


_ 


53w080 


17 21 37.46 


49 55 36.9 


18.22 ±0.05 (8) 


18.37 ±0.09 


17.85 ± 0.05 (8) 


18.00 ±0.09 


_ 


_ 


53w081 


17 21 37.81 


49 57 56.9 


23.99 ±0.13 (1.5) 


23.64 ±0.19 


23.36 ±0.26 (1.5) 


23.01 ±0.31 


_ 


_ 


53w082 


17 21 37.64 


50 08 27.4 


25.01 ± 0.42 (2.5) 


24.86 ± 0.43 


>23.76 


_ 


_ 


_ 


53w083 


17 21 48.93 


50 02 39.8 


22.18 ±0.06 (2.5) 


21.94 ±0.13 


21.52 ±0.06 (2.5) 


21.28 ± 0.13 


_ 


_ 


53w084 


17 21 50.43 


49 48 30.5 


24.78 ± 0.34 (2.5) 


24.61 ± 0.35 


>23.76 


_ 


19.46 ± 0.38 (2.5) 


19.29 ± 0.39 


53w085 


17 21 52.47 


49 54 34.0 


22.17 ±0.06 (2.5) 


22.01 ± 0.10 


21.93 ±0.07 (2.5) 


21.77 ± 0.10 


- 


- 


53w086a 


17 21 56.42 


49 53 39.8 


20.22 ± 0.05 (4) 


20.10 ± 0.08 


19.44 ± 0.05 (4) 


19.32 ± 0.08 






53w086b 


17 21 57.65 


49 53 33.8 


22.08 ± 0.06 (2.5) 


21.69 ±0.12 


20.95 ± 0.06 (2.5) 


20.56 ±0.12 






53w087 






>25.17 




>23.76 




>19.85 




53w088 






>25.17 




>23.76 








53w089 


17 22 01.02 


50 06 51.7 


24.27 ±0.16 (1.5) 


23.84 ± 0.22 


>23.76 




>19.85 




66w009a 


17 18 32.87 


49 55 53.9 


23.11 ± 0.08 (1.5) 


22.68 ± 0.22 


22.36 ±0.11 (1.5) 


21.93 ± 0.24 


16.94 ±0.02 (1.5) 


16.51 ±0.21 


66w009b 


17 18 33.80 


49 56 02.2 


17.71 ± 0.05 (4) 


17.19 ±0.26 


17.19 ±0.05 (4) 


16.67 ± 0.26 


13.79 ± 0.01 (4) 


13.26 ± 0.26 


66w014 


17 18 53.49 


49 52 39.3 


* 












66w027 


17 19 52.11 


50 02 12.7 


18.33 ± 0.05 (8) 


17.99 ±0.19 


17.81 ± 0.05 (8) 


17.46 ± 0.19 






66w031 


17 20 06.87 


49 43 57.0 


22.65 ± 0.07 (2.5) 


22.45 ± 0.12 


22.43 ± 0.10(2.5) 


22.23 ±0.14 


17.96 ±0.10 (2.5) 


17.76 ±0.14 


66w035 


17 20 12.41 


49 57 08.7 


23.47 ±0.11 (2.5) 


23.31 ± 0.14 


23.12 ±0.17 (2.5) 


22.95 ±0.19 


19.10 ±0.29 (2.5) 


18.94 ±0.30 


66w036 


17 20 21.46 


49 46 58.3 


22.79 ± 0.07 (2.5) 


22.60 ± 0.12 


21.82 ±0.07 (2.5) 


21.63 ± 0.12 


17.45 ± 0.06 (2.5) 


17.26 ±0.11 


66w042 


17 20 52.20 


49 42 49.2 


21.21 ± 0.05 (4) 


21.16 ±0.06 


21.01 ± 0.07 (4) 


20.96 ± 0.08 






66w047 


17 21 05.48 


49 56 55.9 


19.30 ± 0.05 (8) 


19.38 ± 0.06 


18.80 ±0.06 (8) 


18.87 ±0.07 






66w049 


17 21 11.21 


49 58 32.9 


22.59 ± 0.07 (2.5) 


22.41 ± 0.11 


22.16 ±0.08 (2.5) 


21.98 ±0.12 






66w058 






>25.17 




>23.76 









the observed galaxies tend to be bluer than the models. This ten- 
dency is not unexpected as the blue rest-frame wavelengths probed 
at those redshifts mean that a smaller amount of recent star forma- 
tion or AGN activity will significantly bluen the galaxy colours. In 
the ir — K) plot, the observed galaxies are in good agreement with 
the models, but again the small number of sources included make 
drawing conclusions difficult. 

The magnitude distribution histograms (Figure are useful 
as they provide a first look at the redshift distribution of the sample 
through the magnitude-redshift relations for radio galaxies (these 
are described in more detail in i]8.U , which indicate that the opti- 
cally faintest objects should lie at the largest distances. The X"-band 
magnitude distribution is ignored here because the low number of 
K-magnitudes taken would not result in a meaningful diagram. 

Both the colour-magnitude and magnitude distribution dia- 



grams are also in reasonable agreement with those shown in Figures 
9 and 10 respectively of Waddington et al. (2000). 



6 SPECTROSCOPIC TNG OBSERVATIONS 

The redshifts for the radio sources identified in the Lynx and Her- 
cules fields are vital in determining their cosmic evolution. Previ- 
ously published spectroscopic or photometric redshifts (Wadding- 
ton et al. 2000 and references therein; Bershady et al. 1994) already 
exist for 19 of the Hercules field sources, and 3 of the Lynx field 
sources had spectroscopic redshifts from the SDSS; the remaining 
sources had no previous redshift information. This section covers 
the spectroscopic observations made of a selection of the sample 
with the multi-object spectrograph DOLORES on the 3.58 m Tele- 
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Table 9. The host galaxy positions and magnitudes for the Lynx field. The radius (in ") of the aperture used for photometry is given in brackets and 3cr limits 
are given for undetected sources. The corresponding radio positions can be found in Table[2] 



Lynx 



Name 


RA 
(J2000) 


DEC 
(J2000) 


r 


r (63.9 kpc) 


i 


i (63.9 kpc) 


K 


K (63.9 kpc) 


SSwl 1 A 

JJW 1 IO 


OS zll ztO 70 
uo 4 J 4U. Iy 


A A 10 9S S 
44 J7 ZJ.J 


99 ni -i- n 07 to s^ 
zz. uj it u.u/ 


9 1 SzL -1-0 19 
Z 1 .04 It U. 1Z 


91 1 f\ 4- 1 (1 S~\ 
Zl.lO IL U.lU (Z.J^ 


90 07 4- n 1 zl 

zu.y / it u. 14 






JJW 1 1 o 


fifi /II Af^ SA 
l/o 4J 40.50 


AA IS AQ 1 
44 J J 4V. / 


91 9Q _l_ n OA (A\ 
L1.L7 It U.UO (4 ) 


91 1A A- 07 

z 1 .Z4 zn u.u / 


90 SO 4- OS (A\ 
ZU.oy IC U.Uo (4) 


90 Szl 4- O OS 
ZU.o4 It U.Uo 






SSwl 90 
JJ W 1ZU 


(\Q A?, S9 87 

UO 4J JZ.O / 


AA 1A 9Q 1 
M-M- LH Ly . 1 


~~~>!A IS 
/ZH.JO 




~^»9^ Af\ 




IS 19 4-0 1(1^9 ^ 


1 7 qa 4- n n 

i / .7\j in u. i j 


JJW 1Z 1 


ns aa nzi m 

Uo 44 U4.U1 


A A *X 1 90 1 
44 J 1 ZU. J 


91 i s _i_ n 1 a (i s 1 ) 
zj.ij in u.i4 \l.j) 


99 no _i_ n \f> 

ZZ. 70 It U. ID 


\91 Af. 
^>ZJ.40 




1 Q is 4- O IS (1 SA 

iy.jj it u.jj \l.j) 


10 1 7 4- n 1A 

iy. i / it u.jo 


cc,„i 99 


ns /i /i 19 in 

Uo 44 1Z. 1U 


A A 11 1 7 S 
44 J I 1 / . J 


9n 7/1 -L- n ns ia\ 

ZU. /4 it U.UJ (4; 


9n fifi -t n n7 
zu.oo it u.u / 


90 s^ 4- n ns (A\ 

ZU.jO it U.Uo y\) 


90 /i 7 4- n no 
zu.4 / it u.uy 






SSwl 91 
JJW 1ZJ 


no AA\A^A 
Uo 44 14.J4 


A A IS no 9 
44 J J UU.Z 


99 on -i- n 1 9 n si 
zz. yu m u.iz (z. j) 


99 71 _|_ n 1 s 
zz. 11 in u. 1 j 


99 QS 4- 1ft (1 s^ 

zz. yo it u.jo V-^--'/ 


99 7Q 4- n 17 
ZZ. 17 It U.J / 


1 7 in 4- O OA (1 SA 

i /.ju zr. u.uo (z.jj 


17 10 4-011 
1/.1U X U.ll 


JJW1Z4 


ns A A 1/1 Q1 
Uo 44 14.yj 


/1/1 IS S9 9 
44 jo JZ.Z 


9 1 99 _i_ n nA //i \ 

Li.LL it U.UO (4) 


91 i< 4. n m 
z 1 . 10 it u.u / 


91 00 4_ n 1 n (A \ 
Zl.Zo It U.lU (4) 


9 1 99 4_ n 1 O 
Z 1 .ZZ It U. 1U 






SSxwl 97 
JJW1Z / 


no AA 97 1 S 
Uo 44 Z / . 1 J 


aa ai ns n 

44 4 J Uo.U 


\a 1 ° -J- n ns f'SA 

14. Io It U.UJ (oj 


1 1 Al _l_ n 90 

IJ.OJ ^ U.Zo 


1 A 1 9 4- O 07 i'S'i 

14. 1Z It U.U/ (0) 


1 1 57 4- n 9S 
1 J.J / It U.Zo 






CC„,1 9Q 

JJWlZo 






\ 9/1 TO 

,>Z4. JO 




>ZJ.40 




\90 1 A 
>*ZU. 10 




SSwl 11 
JJW 1 J 1 


no A A 1 S s 1 
Uo 44 J J. J 1 


AA Af\ V\A 1 
44 40 U4 . 1 


91 1 S -1- 1 S (1 S~l 
ZJ.lo U.1J 


91 09 -I- 1 7 
Z J.UZ x U.l; 


91 70 4- 1 zl /"9 S^ 
Zl. 17 It U. 14 (Z.JJ 


91 fO 4- O 1 f\ 
Z 1 .OZ It u. 10 






CC„,1 "39 






\9/1 IS 
^>Z4. Jo 




-^91 Afx 
>ZJ.40 




\90 1 A 
>ZU. 10 




JJWl JJ 


nS AA 17 9/1 
Uo 44 J / .Z4 


a a 9 A nn /I 

44 ZO UU.4 


9S S1 -t 1 01 / 1 S 1 ! 
ZJ.Jl it l.UJ \l. J) 


9S K J. 1 OS 
ZJ. 1 J it l.UJ 


\91 /1A 
,>ZJ.40 




\in no 

s> iy.ys 




jjWIjj 


ns A A a\ in 

Uo 44 4 1 . 1 U 


AA 91 17 7 
44 Z 1 j 1 . 1 


>Z4. Jo 




\91 A£, 
>Zj.40 




1 1 91 4_ n n9 ( i n\ 

Ij.Zj it U.UZ (1UJ 


1 1 1 9 4_ O 91 
1 J. 1Z it U.ZJ 


SSwl If. 

jjWIjO 


ns a A a s no 

Uo 44 4J.UV 


AA 19 97 1 
44 JZ LI .J 


91 sn 4- n 99 ( 1 s^ 

ZJ.oU it U.ZZ (l.jj 


91 A^ 4- O 9S 
ZJ.4J It U.Zo 


\91 A£. 
>ZJ.40 




1Q 1 7 4- O 1 A / 1 S\ 

iy.i / itu.io^i.j^ 


1 S S 1 4- O 9/1 
lo.ol it U.Z4 


JJWl J 1 






\9/1 IS 
^>Z4. Jo 




\91 /1A 
^•Zj.40 








SSwl IS 
JJW 1 Jo 


nS AA Szl AS 
Uo 44 J4.4J 


zlzL 9A 99 n 
44 ZO ZZ.U 


\ 9/1 10 




\91 Af. 
^>ZJ.40 




1Q 79 4- n 9S (\ SA 
Y7.IL It U.ZJ (l.J^ 


10 1/14-0 11 

iy.j4 it u.ji 


JJW14U 


ns /is nA nA 

Uo 4j UO.UO 


/i/i /in /i i 9 

44 4U 41 .Z 


9n 79 -U- n ns //i\ 

ZU. /Z it U.Uj (^4^ 


90 7S 4- O OS 
ZU. / J it U.Uj 


90 Q/i 4- n ns 

ZU. y4 it U.Uo (4) 


90 4- n ns 
zu.yo it U.Uo 






jjW141 






\9 ,1 TO 

_>Z4. Jo 




\91 /1A 
>Zj.40 








J jW 14 Ju 


ns /is ns a 9 

Uo 4 J UJ.OZ 


/1/1 9S /19 Q 
44 ZJ 4Z.y 


9^ is -t n Q9 <"1 s^ 
zj.jo it u.yz (i.jj 


9^ ni 4- 1 on 

Zj.Uj it l.UU 


■^91 Afx 
> ZJ.40 




s> ly.oj 




JJW14JD 


ns /IS CiA 9S 
Uo 4j U4.ZJ 


AA 9S SI 1 
44 ZJ J J. J 


9^ A£, -i- n OS (\ S~i 

zj.40 it u.ys (i.jj 


95 1 1 4. HQS 

zj. 1 1 it u.yo 


•\91 Afx 
> ZJ.40 




>■ ly.oj 




sSw 1 An 

JJWI4/ 


ns /is 91 si 

Uo 4j Zj.oJ 


AA SO 9/1 A 
44 jU Z4.0 


91 n7 4. n n fi ^\ 

Zj.U/ it U.lJ (Z.J) 


99 on 4_ n 1A 
zz. yu it u. 10 


-^91 AA 
>Zj.40 




1 7 AS 4- n n7 /9 s\ 
1 / .Do it U.U/ (Z.Jj 


1 7 s i x n 1 1 

1 / .jl it U. 1 1 


^ s^t i /io 


ns a c; T7 i 7 
Uo 4 J z / . 1 / 


/1/1 ^5 nc n 

44 J J zj.y 


1 ^ aq -L n ns /q\ 
10.4V it U.UJ (0) 


1 a ia 4_ n no 
10. 34 it u.uy 


1 ^ Q/i 4_ n m /'s^ 
lJ.y4 it U.U/ (0) 


1 ^ s n 4_ n 1 n 
1 j.oU it U.lU 






JJWl JU 


nS /IS 9Q /1 7 

Uo 4j zy.4 / 


aa so 17 /I 
44 JU J / .4 


9n S9 -t n ns //i\ 

ZU.oZ it U.Uj (4; 


9n 7n 4- n ns 

ZU. /U it U.Uo 


9n 1 9 4- n n7 

ZU.1Z it U.U/ (4) 


on nn 4- n no 
zu.uu it u.uy 






JjW 1 j4 


ns a s a i in 

Uo 4j 4 1 . JU 


AA /in 1 1 Q 
44 4U 1 1 .y 


10 io 4- n ns /s\ 
iy. 17 it U.Uj yo) 


1 o 9 ^ 4_ n n^ 
ly.zj it u.uo 


1 s so 4- n ns <"S~i 
lo.jy it u.uo \<s) 


1 s A/i 4- n ns 

10.O4 it U.Uo 






JJWl JJ 






_>Z4. Jo 




\91 A£, 
>ZJ.40 




\1Q o c 

s> ly.oj 




JjWIjO 


ns /is sn Q9 

Uo 4j jU.VZ 


AA 1Q S 1 S 
44 37 J 1 .J 


99 7^ 4_ n in 1*9 s^ 
ZZ. / J it U.lU (Z.J J 


99 s^ 4- n 1 /i 

ZZ. JO it U. 14 


91 ni 4- O 17 /9 S\ 
Zj.Uj it U.j/ \L.j) 


99 S/1 4- n 18 
ZZ. 54 it U.Jo 


17 97 4_ n ns /9 s\ 

1 1 .LI it U.UJ (Z. j ) 


1 7 n7 x n 1 1 

1 /.U / It U. 1 1 


55wl57 


08 46 04.44 


44 45 52.7 


22.07 ± 0.07 (2.5) 


21.77 ±0.16 


21.45 ± 0.11 (2.5) 


21.15 ±0.18 






55wl59a 


08 46 06.67 


44 51 27.5 


23.54 ± 0.08 (2.5) 


23.38 ± 0.22 


>23.46 








55wl59b 


08 46 06.66 


44 50 53.8 


18.70 ±0.05 (8) 


18.74 ib 0.05 


18.11 ± 0.07 (8) 


18.15 ±0.07 






55wl60 


08 46 08.57 


44 36 47.4 


21.40 ±0.06 (4) 


21.33 ± 0.07 


20.24 ± 0.07 (4) 


20.17 ±0.08 






55wl61 


08 46 27.46 


44 29 57.1 


20.07 ± 0.05 (4) 


19.94 ±0.08 


19.46 ± 0.07 (4) 


19.33 ±0.10 






55wl65a 


08 46 34.78 


44 41 37.6 


21.36 ±0.06 (4) 


21.31 ±0.06 


20.31 ± 0.07 (4) 


20.26 ± 0.07 






55wl65b 


08 46 33.37 


44 41 24.4 


21.65 ±0.07 (4) 


21.62 ±0.07 


20.89 ± 0.08 (4) 


20.86 ± 0.08 






55wl66 
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44 30 53.5 


22.72 ±0.10 (2.5) 


22.54 ±0.13 


21.99 ±0.16 (2.5) 


21.81 ± 0.18 
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20.78 ± 0.08 
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>24.38 




>23.46 








60w039 


08 44 42.50 


44 45 32.5 


17.23 ± 0.05 (8) 


17.08 ± 0.08 


16.78 ± 0.07 (8) 


16.63 ±0.10 
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08 45 14.00 


44 53 08.7 


21.85 ±0.06(2.5) 


21.63 ± 0.12 


20.79 ± 0.08 (2.5) 


20.57 ±0.13 






60w067 






>24.38 




>23.46 








60w071 


08 46 00.34 


44 43 22.1 


23.44 ±0.18 (2.5) 


23.28 ± 0.20 


>23.46 








60w084 


08 46 40.23 


44 33 44.7 


17.79 ± 0.05 (8) 


17.58 ±0.11 


17.06 ± 0.07 (8) 


16.85 ± 0.12 







scopio Nazionale Galileo (TNG), along with the redshift estimation 
methods used for the remaining sources. 

DOLORES, the Device Optimized for the LOw RESolution, 
consists of one Loral back-illuminated and thinned 2048x2048 
pixel CCD, with a scale of 0.275"/pixel, resulting in a field of 
view of 9.4'x 9.4'. For multi-object spectroscopy (MOS) observa- 
tions, rectangular masks with dimensions 6.0'x 7.7', are used. Ver- 
tical slits of constant width (either 1.1" or 1.6") but varying length 
are drilled in the masks according to the positions of the various 
sources for which spectra are required. 

Since some uncertainties in the astrometry remained the 1.6" 
slit was used for the observations. Only 10 masks were permitted 
per observing run, so it was decided to create 5 masks of varying 
position angle for each for the two fields, thus covering as many 



sources as possible. £]6. II below describes the mask creation and 
source selection process. 

The DOLORES observations took place from 18th to 20th 
April 2004. All nights were photometric and the standard star Feige 
34 was observed (with a 5" longslit) at regular intervals. Each mask 
observation consisted of several long exposures to avoid satura- 
tion of the CCD, and to allow cosmic ray hits to be identified in 
the final spectra. The good weather conditions also allowed two 
observations with the 1.5" longslit of three sources (two in Her- 
cules and one in Lynx), which were not included on the masks. All 
observations were carried out using the LR-R grism which has a 
wavelength range of 4470-10360A and a resolution of 11. OA. Full 
details of the observations can found in Table 1111 

The masks for the Lynx field were, on average, observed for 
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Figure 3. Colour magnitude diagrams and completeness limits for the Lynx (open diamonds; dot-dashed lines) and Hercules (rilled diamonds; dashed 
lines) results. Overplotted is the predicted behaviour of an L* type galaxy for two models (Bruzual & Chariot, 2003): passive evolution following an 
instantaneous burst of star formation at z = 5 (dotted line), and exponentially declining star formation, formed at z = 10, with an e-folding time of 1 
Gyr (dot-dot-dot-dashed lines). 
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Figure 4. Magnitude distribution for the combined Lynx and Hercules complete sample in / (left) and i' (right). The shaded regions show the undetected 
sources at the magnitude limits. 



less time than those for Hercules due to the early setting time of the 
Lynx field. 



6.1 MOS Mask Creation 

The mask limitations meant that not all the radio sources could be 
included in the observations, therefore a ranking system was in- 
troduced to ensure that the optimum number of interesting, high- 
redshift, sources were included. First, r-band magnitudes were es- 
timated for all the sources using a rough m zpt and ignoring the 
airmass and extinction corrections. (These magnitudes could only 
be estimated as, at the time these observations were being prepared, 
only the non-photometric April 2003 WFC data were available.) 

The resulting magnitudes were then used to give an indication 
of the redshifts of the sources through the r-z relation (Snellen et al. 
1996). Those sources with previously published redshifts were ob- 
viously excluded from the rankings. Table 1 101 gives the full details 
of the ranking scheme. 

The masks were then arranged such that the maximum num- 
ber of first rank sources would be observed. This was done using 
an IDL script which allowed possible combinations of objects to 



Table 10. The ranking scheme for the radio sources. 



Redshift 


r magnitude 


Rank 


z>0.5 


21-22 


1 




>22 


2 


z<0.5 


<21 


3 




<15 


4 



be rotated to determine the best fit parameters for each mask. Table 
[TT]gives the final determined parameters. In the Lynx field 56% of 
all the sources were included, and of these 30% were first rank and 
78% were second rank or above. The Hercules field was slightly 
better with 62% of the total number included; 80% classed as sec- 
ond rank or higher and 11% first rank. 



6.2 Data reduction 

The initial data reduction of the DOLORES observations was very 
similar to the WFC reduction process; again the IRAF software 
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Table 11. Summary of the TNG observations. H and L refer to the Hercules and Lynx fields respectively. 



Mask 


PrwiHnn Anolp (°\ 

l VJfilLlUll ^\.lltlt^ ^ } 


Centre (720001 


Slit I enoth ("1 


Ohcprvnti mi 

V-/Uc>d V dll Wll 




Seeinp ("1 
seeing v ) 










Date 


Time 




Ll 





8 45 53.78+44 43 01.8 


12 


19/04/04 


4x1 500s 


1.3 


L2 


20 


8 45 24.20 +44 53 43.4 


12 


20/04/04 


3x1 800s 


0.8 


L3 





8 44 55.98 +44 44 04.8 


12 


19/04/04 


2x1 800s 


1.3 


L4 


-5 


8 44 36.14+44 46 50.1 


10 


18/04/04 


3x1200s 


1.3 


L5 


-95 


8 44 00.00 +44 37 08.0 


12 


18/04/04 


3x1500s 


1.3 


HI 


-40 


17 21 51.57+49 50 31.3 


11 


18/04/04 


4x1 500s 


1.3 


H2 


10 


17 21 52.74 +50 09 11.4 


12 


19/04/04 


4x1 800s 


1.3 


H3 


-15 


17 20 14.86 +49 47 20.3 


12 


18/04/04 


3x1800s, 1x492s 


1.3 


H4 


-45 


17 20 09.00 +50 01 33.6 


12 


20/04/04 


4x1 800s 


0.8 


H5 


-90 


17 19 05.97 +49 45 59.4 


12 


19/04/04 


3x1800s, 1x100s, 1x1200s 


1.3 



package was used throughout. First a master bias was constructed 
for each night by averaging all the appropriate bias frames together; 
this was then subtracted from the rest of the exposures. The science, 
flat-field and arc data were then grouped according to mask. The 
calibration arc-lamp used was HeNeAr. 

The science and arc images for each mask were then combined 
and split into individual two-dimensional spectra for each source. 
The normalised flat field for each source was then applied and the 
strong background sky-lines were removed. This process was not 
perfect and some background residuals remained; these had to be 
taken into account when attempting to identify spectral features. 

One-dimensional spectra were then extracted for the sources 
and their corresponding arcs. Lines in the arc spectrum were then 
identified and the resulting calibration was applied to the science 
spectra. The calibration was improved by adding sky-lines ([OI] 
at 6300.3A and 5579A and Na at 5896A), extracted from the un- 
background subtracted science spectra, to the arc spectra to im- 
prove the wavelength coverage. 

The spectra were then flux calibrated using observations of a 
standard star, Feige 34. The standard star spectrum was also used 
to correct, to a first approximation, for sky absorption features. 

6.3 Spectroscopic Results 

The spectroscopic observations included 41 sources in total; 17 in 
the Hercules field and 24 in the Lynx field. The resulting spec- 
tra yielded 3 and 1 1 definite redshifts in the Hercules and Lynx 
fields respectively. A single line was detected in a further 4 spectra 
in Lynx and 2 in Hercules; this mainly provided a redshift 'best- 
guess' only, except where the line identification was obvious (e.g. 
broad Mgll). 65% of the redshifts were from emission lines and 
the remainder were from absorption lines only. These results are 
given more fully in Table [13] and the spectra for sources with one 
or more detected lines or absorption features, are shown in Figure 
|C1| The list of sources observed with DOLORES, but for which no 
lines were detected can be found in Table [12] Notes on individual 
sources can be found in Appendix I A2 1 



7 CLEANING THE SAMPLE: IDENTIFYING QUASARS 
AND STARBURST GALAXIES 

Not all the radio sources detected in the two fields will be radio 
galaxies; some contamination of the sample by quasars and star- 
burst galaxies is inevitable. The quasars need to be identified as 



Table 12. The sources targeted in the DOLORES observations with no lines 
detected. 



Mask 


Source 


HI 


53w084, 53w086a, 66w058 


H2 


53w082, 53w087, 53w088, 53w089 


H3 


53w069 


H4 


66w035 


H5 


53w054a, 53w057, 53w061 


Ll 


55wl56, 60w071 


L2 


55wl47 


L3 


55wl38 


L4 


55wl27, 55wl32, 60w032 


L5 


55wll8, 55wl23 



the photometric redshift estimation methods are not valid for them, 
whilst the starburst galaxies need to identified and removed. 

7.1 Identifying starburst galaxies 

The radio emission from starburst galaxies is mainly due to syn- 
chrotron emission from supemovae instead of, as in radio galax- 
ies, accretion onto a supermassive black hole. The radio luminos- 
ity function of Best et al. (2005) shows that, in general, the radio 
power of these galaxies is lower than that of AGN and that it is only 
at these very low radio powers (~ 10 23 WHz~ 1 ) that their number 
density dominates; this suggests that only low-power sources need 
to be considered here. There were 16 radio sources in the sample 
which have P rad «S 10 24 WHz" 1 and were therefore possible 
starburst galaxy candidates. 

The main way to distinguish between radio and starburst 
galaxies is through examination of the emission line ratios of their 
spectra. Kauffmann et al. (2003) classify a source as an AGN if 

> log ( /[NII] /° / (L))-o.o5 + L3 ' (3) 

where /[OIII], /(H/3), /[Nil] and /(Ha) are the fluxes of the 
respective emission lines; [OIII] (5007A), [Nil] (6583A), Hp 
(486lA) and H a (6563A). This classification method obviously 
requires a source to have spectroscopic data with the right lines 
detected. Of the 16 candidates, 11 have DOLORES spectra but 
none have all four of the necessary emission lines. However, 3 can- 
didates (66w027, 55wl37 and 60w039) have both [Nil] and Ha 
detected and one, 55wl50, has both [OIII] and H/3, hence an in- 
dication of their classifications is possible; the results of this are 
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Table 13. Spectroscopic redshifts and line information for the Hercules and Lynx fields. 



Mask 


Source 


A 

(A) 


Line 


Flux 

(x 10~ 16 crg s _1 cm -2 ) 


^fwhm 

(kms -1 ) 


W 

(A) 


z 


Final z 


Hercules 


H4 


53w070 


6480.5 


Mgll 


0.48 ± 0.07 




27 ±5 


1.315 ± 0.001 


1.315 ±0.001 


HI 


53w086b 




4000A break 








0.73 ± 0.01 


0.73 ± 0.05 


H4 


66w027 


7136.6 


Ha 


16.9 ± 2.55 


1091 ± 228 


39 ±4 


0.087 ± 0.001 


0.086 ± 0.002 






7136.5 


Nil 


6.9 ± 2.55 


1074 ± 228 


38 ±4 


0.084 ± 0.001 








7311.7 


[sm 


~3.94 




~6 


0.088 ± 0.001 








6401 


NaD 
















5282 


H/3 


10.77 ± 1.94 




17 ±3 






H3 


66w031 


6751.8 


[Oil] 


2.47 ± 0.25 


987 ± 251 


508 ± 205 


0.812 ±0.001 


0.812 ± 0.001 






8822 


H/3 












H3 


66w036 


8273.4 


G-band 








0.924 ± 0.001 


0.924 ± 0.001 






7581.4 


CaH 








0.927 ± 0.012 








7650 


CaK 












Lynx 


L5 


55wll6 


7281.5 


CaH 


_ 


_ 


_ 


0.854 ±0.031 


0.851 ±0.007 






7343.4 


CaK 


_ 


_ 


_ 


0.851 ±0.008 




L5 


55wl24 


6536.5 


Mgll 


4.85 ± 0.49 


5471 ±894 


39 ±4 


1.335 ±0.003 


1.335 ±0.003 


L4 


55wl28 


8159.3 


[Oil] 


0.17 ±0.03 


_ 


_ 


1.189 ± 0.001 


1.189 ±0.001 


L4 


55wl31 


7917.1 


[Oil] 


0.75 ±0.10 


409 ± 281 


_ 


1.124 ±0.001 


1.124 ±0.001 


L3 


55wl37 


5755.8 


[OIII] 


1.64 ±0.22 


_ 


4± 1 


0.150 ±0.001 


0.151 ±0.001 






5703.8 


[OIII] 


9.94 ± 2.56 


_ 


26 ±7 










7567.1 


Ha 


11.7 ±2.63 


2000 ± 266 


37 ±4 


0.153 ± 0.001 








7575.2 


[Nil] 


13.1 ±2.56 


1975 ± 265 


36 ±4 


0.150 ±0.001 








7738.4 


[sm 


2.77 ±0.31 


- 


- 


0.152 ±0.001 








6778.4 


NaD 


_ 


_ 


_ 


0.150 ±0.001 








7248.8 


[OI] 


0.57 ±0.15 


_ 


_ 


0.151 ± 0.001 




L3 


55wl40 


7514.3 


Mgll 


12.10 ± 1.40 


6192 ± 2647 


36 ±4 


1.685 ±0.012 


1.685 ±0.012 


L2 


55wl49 


6793.7 


NaD 


_ 


_ 


_ 


0.152 ±0.001 


0.151 ±0.001 






5950.6 


Mgb 


- 


- 


- 


0.150 ±0.001 








5747.7 


[OIII] 
















5601.2 


H/3 












L2 


55wl50 


7359.3 


[OIII] 


6.16 ±0.64 


727 ± 233 


67 ± 8 


0.470 ± 0.001 


0.470 ± 0.001 






7292.6 


[OIII] 


1.65 ± 0.21 




22 ±3 


0.471 ± 0.001 








7138.4 


H/3 


0.96 ±0.16 




21 ±4 


0.469 ± 0.001 








9660.4 


Ha 












LI 


55wl54 


6475.6 


H/3 








0.332 ± 0.001 


0.330 ± 0.001 






6891.1 


Mgb 








0.332 ± 0.001 








5717.9 


G-band 








0.330 ± 0.002 








5227.3 


CaH 








0.329 ± 0.001 








5277.5 


CaK 
















5456.5 


B8 












LI 


55wl57 


6760.2 


H7 








0.558 ± 0.001 


0.559 ± 0.002 






6201.4 


CaK 








0.563 ± 0.001 




LS 


55wl60 


6292.3 


CaH 








0.600 ± 0.002 


0.600 ± 0.002 






6352.0 


CaK 












L5 


60w016 


7237.0 


CaH 








0.840 ±0.001 


0.840 ±0.001 






7332.7 


CaK 












L5 


60w024 


6609.9 


[Oil] 


0.41 ± 0.06 


814 ± 374 


13 ±2 


0.774 ± 0.001 


0.774 ± 0.001 






6974.2 


CaH 
















7032.2 


CaK 












L3/4 


60w039 


7562.2 


Ha 


51.0 ±6.18 


939 ± 223 


47 ± 5 


0.152 ±0.001 


0.151 ±0.001 






7585.2 


[Nil] 


11.4 ±6.12 


925 ± 223 


46 ±5 


0.149 ±0.001 








7738.7 


[sm 


8.32 ± 1.33 






0.152 ±0.002 








6786.8 


NaD 








0.151 ± 0.001 








5597.6 


H/3 








0.152 ±0.001 




L2 


60w055 


6405.0 


[Oil] 


1.40 ±0.14 


529 ± 255 


45 ±5 


0.718 ±0.001 


0.718 ±0.005 






6754.4 


CaH 








0.717 ± 0.001 








6818.8 


CaK 
















7397.9 


G-band 
















7437.8 


H7 
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outlined below. The fluxes for these lines can be found in Table 
n~3l A further two candidates, 55wl35 and 60w084, were included 
in the spectroscopic observations of the SDSS. The resulting spec- 
trum for 55wl35 clearly shows it to be a starburst galaxy since 
/(Ha) > /[Nil] and /(H/3) > /[OIII] whereas the spectrum 
for 60w084 suggests that it is an AGN since /(Ha) ~ /[Nil] and 
/[OIII] » /(H/3). 

• 66w027 (z=0.086, Pi. 4 ghz=10 22 00 WHz- 1 ) 
log(/[NII]//(Ha)) = —0.38 for this source which suggests that 
it is a starburst galaxy but is not sufficient to be unambiguous. 
However, the H/3 line was also detected whilst the [OIII] line was 
not, implying that log(/[OIII]//(H/3)) < 0.0 and that this source 
is a starburst galaxy. 

• 55wl37 (z=0.151, Pi.4GHz=10 22 ' 97 WHz- 1 ) 
log(/[NII]//(Ha)) = 0.05 for this source which places it 
firmly in the AGN region. This is further supported by the 
indication of extended radio emission visible in the radio image. 

• 60w039 (z=0.151, Pi.4GHz=10 22 - 58 WHz- 1 ) 
log(/[NII]//(Ha)) = —0.65 for this source which places 
it firmly in the starburst region. 

• 55wl50 (z=0.470, Pi.4GHz= 
log(/[OIII]//(H/3)) = 0.78 for this source which, coupled 
with the high radio power, strongly suggests that it is an AGN. 



zlO^^WHz -1 ) 



100.00 r 



10.00 r 



1.00 : 



0.10 



^0w0Z4 



<£0v»055 



55w160 

* 55w15Z- 



1 00 



1 000 

.4CHz flux density / /xJy 



10000 



Figure 5. Radio-[OII] flux relationships for starburst galaxies (solid line) 
and radio galaxies (dashed line). The dotted line shows the effect of two 
magnitudes of dust extinction in the measurement of the star formation rate. 
Overplotted are the starburst candidate sources for classification. 



The five remaining starburst candidates with spectra can be 
classified by comparing their /([Oil]) and radio fluxes. Best et al. 
(2002) derive a rough relationship between these two fluxes for 
starburst galaxies by equating the theoretical correlation between 
mean star-formation rate and [Oil] luminosity (Barbaro & Pog- 
gianti, 1997) with that for radio (Condon & Yin, 1990) giving 



Si. 4 



fxJy 



11.0 



/[OH] 



10 _16 ergs _1 cm~ 



(4) 



where iSi.4GHz is the radio flux of the object in question at a fre- 
quency of 1.4 GHz and z is the redshift. It should be noted that dust 
extinction can cause the measurement of the star-formation, from 
the [OIII] flux, to be underestimated. 

The radio-[OII] flux relationship for AGN is an extrapolation 
of the results of Willott et al. (1999) for powerful radio galaxies to 
the sub-mjy levels of these sources. This gives, after converting to 
1.4 GHz (assuming a spectral index of 0.8), 



1.4GHz 

fijy 



■■ 4.7 x 10 4 



/[OH] 



. 10 _16 ergs _1 cm~ 



(5) 



It should be noted however, that AGN show considerable scatter 
around this line. 

Equations [4] and [5] are shown in Figure [5] Overplotted are the 
[Oil] and radio fluxes for 60w055 and 60w024. The remaining 3 
sources (60w016, 55wl60 and 55wl57) lack an [Oil] detection so 
an upper limit for their [Oil] flux is plotted instead. It is clear from 
inspection of this diagram that all these sources lie well below the 
starburst region, suggesting that they are AGN. This is in accor- 
dance with their radio powers, all of which are > lO^^WHz" 1 . 

This leaves 5 candidates still to be classified. Two (66w009b 
and 55wl27) have radio powers that are < 10 23 ' 01 WHz _1 which 
strongly suggests that these are also starburst galaxies. The final 
three candidate sources (55wll8, 55wl22 and 55wl61) all have 
P > 10 23 WHz _1 which is comparable with the powers of candi- 
dates already identified as AGN. Therefore, on the balance of prob- 
ability, these three objects are classified as AGN also. 

In summary, therefore, there are five starburst galaxies which 



need to be removed from the sample: 66w027, 66w009b, 55wl27, 
55wl35 and 60w039. 



7.2 Identifying quasars 

Possible quasars in the sample were identified in three ways; firstly 
on the basis of their point-like appearance (determined by mea- 
suring the FWHM using gaia) in the r-band images, secondly by 
looking for broad lines (3>1000kms _1 ) in the available spectra and 
thirdly by examining their optical (r — i) colours. The quasar can- 
didates selected in this way are described below. 

• 53w061 - This source has a pointlike appearance (FWHM of 
1.53"in r-band, c.f. seeing of 1.5") and a blue (r — i) colour of 
0.36. It is also classed as a likely quasar (Q?) by Rron et al. (1985) 
on the basis of its colour. It was included in the DOLORES obser- 
vations but no continuum was detected. It is therefore classed as a 
quasar here. 

• 53w065 - This source has a blue (r — i) colour of -0.31 but 
Waddington et al. (2001) classify it as a galaxy on the basis of the 
narrow lines seen in its spectrum. 

• 53w075 - This source is classified as a quasar by Rron et al. 
(1985) on the basis of its spectrum; this is supported by its pointlike 
appearance. 

• 53w080 - This source is classified as a quasar by Rron et al. 
(1985) on the basis of its spectrum; this is supported by its pointlike 
appearance. 

• 55wl21 - This source is classified as a quasar by Rron et al. 
(1985) on the basis of its colour. However, whilst its non-detection 
in i does suggest a blue colour, this is not supported by its, not very 
blue, (r — K) of 3.8, and it does not appear to be pointlike. It is 
therefore classed as a galaxy here. 

• 55wl24 - This source is classified as a quasar on the basis of 
its very broad Mgll line and blue (r — i) colour of -0.06. 

• 55wl40 - This source is classified as a quasar on the basis of 
its very broad Mgll line and very blue (r — i) colour of -0.21. 
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In summary, therefore, the objects classed as quasars in the 
sample are 53w061, 53w075, 53w080, 55wl24 and 55wl40. 



8 REDSHIFT ESTIMATION AND RESULTS 

Redshifts have now been determined for 44% of sources in the 
Lynx field and 62% of sources in the Hercules field, either through 
the observations described above or with previously published re- 
sults. For the remainder, estimated redshifts were calculated instead 
using the two different magnitude-redshift relations, K-z and r-z, 
outlined in ij8. II below. The K-z relation, the more accurate of the 
two, was used for all sources detected in the UKIRT observations; 
the remaining sources were estimated with the r-z. 

















1 ■ 1 






Photomelri 


'z - ZsKmlll z) vs z 



8.1 The K-z and r-z magnitude-redshift relationships 

The K-z relation is a tight correlation between the A"-band magni- 
tudes and redshifts of radio source host galaxies; it exists because 
the radio host galaxy population is made of up of passively evolv- 
ing, massive elliptical galaxies. 

At high redshift the relation is slightly different for the radio 
surveys of different flux density limits (Willott et al. 2003), which 
therefore implies that it is dependent on the flux density of the radio 
source. The relation used here is that found for the 7C survey as 
the source fluxes for that sample are the best comparison to those 
considered here. The relation itself is given by Willott et al. (2003) 
as, 

K = 17.37 + 4.53 log z - 0.31(log zf . (6) 

The majority of the sources with a host galaxy detection but 
without redshifts, were also not included in the UKIRT observa- 
tions so, therefore, the K-z relation outlined above cannot be used 
for them. Instead, since all these sources have an r-band magni- 
tude, the r-z relation was used. This relation describes a correlation 
between the redshift and r-band magnitude of the host galaxies 
of Gigahertz Peaked Spectrum (GPS) radio sources (Snellen et al., 
1996). 

In general, redshifts estimated using the r-z relation are not as 
reliable as those found using the K-z, especially above z ~ 0.6 
where the r-band samples shortward of the 4000A break. At these 
bluer wavelengths, factors such as star formation or, in more pow- 
erful radio galaxies, the 'alignment effect' (where the optical struc- 
tures of the galaxy align with the radio jets (e.g. McCarthy et al., 
1987)) can result in considerable scatter in the measured magni- 
tudes. GPS sources are less affected by this scatter (c.f. Figure 6 
of Snellen et al. 1996) so are ideal for the redshift estimation here. 
The GPS relation is 

r g = 22.7 + 7.4 log z. (7) 

and is valid up to z ~ 1; at higher redshifts it is unreliable due to 
the scarcity of measurements and the blue rest-frame wavelength 
range it is sampling. 

The relation given in Equation [7] above assumes that the r— 
band magnitude used was observed using the Thuan-Gunn filter 
system. It therefore had to be transformed to the Sloan filter sys- 
tem, via the Johnson filter system, before it could be used. In the 
following description capital letters indicate the Johnson system, 
lower-case letters indicate the Sloan system and g subscripts indi- 
cate the Thuan-Gunn system. To convert from Johnson to Thuan- 
Gunn, Jorgensen et al. (1994) give: 



Figure 6. A comparison of the different methods used to obtain redshifts. 
Filled diamonds indicate the comparison between the two redshift estima- 
tion methods whilst empty diamonds and triangles indicate the comparison 
between the spectroscopically determined redshifts and the r-z and K-z 
estimates respectively. Square outlines indicate redshifts that came from 
single-line spectra. The dotted lines are the \a values of the spectroscopic 
Az, ±0.15. 

r g = .R + 0.111( 5g -r g )+ 0.317, (8) 

whilst to convert from Johnson to Sloan, EquationQjfrom Smith et 
al. (2002) gives: 

r' = R + 0.16{V - R) + 0.13. (9) 

Making the assumption that (V — R) ~ (g s — r g ) ~ 1 for radio 
galaxy hosts, then r = r g — 0.14 and the r-z relation becomes 

r = 22.56 + 7.4 log z. (10) 

The magnitude-redshift relations need aperture corrected 
magnitudes to give reasonable results, but the estimated redshift is 
needed to perform the aperture correction (as described in i]4.3.4| l. 
To solve this problem an IDL script was written to iterate the red- 
shift estimates and subsequent aperture corrections until they con- 
verged on a final value. 

Table [14] summarises all the redshift information for the 
sources in the two fields. The quoted redshift for each source is ide- 
ally a spectroscopic or previously published one, if neither of those 
is available then the redshift is one from a K-z estimation and fi- 
nally, for sources with no A"-magnitude or spectrum, the redshift 
given is an r-z estimate. 

8.2 Redshift comparison 

A valuable test of the redshift estimation comes from comparing the 
estimates with the DOLORES MOS and the previously published 
(Waddington et al. 2000; 2001) spectroscopic and photometric red- 
shifts. The results of this comparison for both r-z and K-z redshift 
estimates are shown in Figure[6j in general the agreement between 
the spectroscopic and estimated results is very good with a la Az 
of 0.15. The agreement is also reasonably good for the photometric 
redshifts. 

Additionally, several of the sources in the sample have both 
r and K-band magnitudes, thus providing a useful means of com- 
paring the two methods of redshift estimation. The solid diamonds 
plotted in Figure|6]show the difference in the two estimates ((K-z)- 
(r-z)) plotted against the K-z value. The two relations give similar 
redshifts up to Zk z ~1.5, but, for the redshifts higher than this, the 
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Table 14. The redshifts found for the sources in the complete sample, 
(la) and (lb) indicates a previously published value, (a - spectroscopic 
(Waddington et al. 2000 and references therein; Bershady et al. 1994), b— 
photometric (Waddington et al. 2001)), (lc) indicates the redshift came 
from the SDSS, (2) indicates a DOLORES spectroscopic value, (3) is a 
K-z estimate, (4) is a r-z estimate. 



Hercules 






Lynx 






Name 


z 


Origin 


Name 


z 


Origin 


53w052 


0.46 


la 


55wl 16 


0.851 


2 


53w054a 

J J W UJT'u 


1.51 


3 


55wl 1 8 

JJW 1 1 C5 


0.66 


4 


53w054b 


3.50 


3 


55wl20 


1.35 


3 


53w057 

JJ W v/J / 


1.85 


4 


55wl21 


2.57 


3 


53w059 


1.65 


4 


55wl22 


0.55 


4 


53w061 


2.88 


lb 


55w123 

J J W 1 J 


0.87 


3 


53w0fi? 


0.61 


la 


55wl24 


1.335 


2 


53w065 

JJW \J\J*J 


1.185 


la 


55wl27 


0.06 


4 


53w066 


1.82 


lb 


55wl28 


1.189 


2 


53w067 

JJW \JVJ 1 


0.759 


la 


55wl31 


1.124 


2 


53w069 

JJW \J\Jy 


1.432 


la 


55wl32 


>4.4 


3 


53w070 

J J W \J 1 \J 


1.315 


2 


55wl 33 

JJW 1 JJ 


2.24 


4 


53w075 

J J W \J 1 J 


2.150 


la 


55wl 35 

J J W 1 J J 


0.090 


lc 


53w076 

J J W \7 / u 


0.390 


la 


55w1 36 

JJW 1 JV.J 


2.12 


3 


53w077 

J J W \^ / / 


0.80 


la 


55wl37 


0.151 


2 
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K-z value is much greater than the r-z further suggesting the lower 
accuracy of the r-z relation at these values. 

8.3 The redshift distribution of the sample 

Now that redshift information has been obtained or estimated for a 
large proportion of the objects in the two fields, redshift histograms 
can be constructed to compare the radio source distributions for the 
two fields; these are shown, split into the different redshift methods, 
in Figures [7(a)land [7(b)1 Both histo grams, peak at redshifts before 

2=1.0. 

These histograms also illustrate the lack of definite redshifts 



Hercules 




1 2 3 4 5 



(b) 

Figure 7. The redshift distribution for sources in the Hercules (a) and Lynx 
(b) fields subdivided to show the contribution from the different redshift 
methods. 

for the Lynx field, especially at the high end, compared to the Her- 
cules field. Whilst many of the Lynx sources were included in the 
DOLORES MOS observations, lines in the resulting spectra tended 
to be detected in the brighter, lower redshift objects. The high-z 
end of both fields, however, is populated mainly by sources with 
less accurate redshift estimates. 



9 CONCLUSIONS 

In summary, the data presented here form a complete sample of 
81 radio sources above the limiting flux density of Si.4GHz > 0.5 
mjy. Of these, only 12 remain unidentified after the optical and 
infra-red observations. Redshifts either previously existed or have 
been determined for 49% of the sample; the remaining objects with 
a host galaxy detection, have redshift estimates. 

The radio, optical, infra-red and MOS observations outlined 
here define a multi-wavelength sample which will form the ba- 
sis for investigations into the high-redshift space density and sub- 
sequent cosmic evolution, if any, of FRI radio galaxies. Paper II 
will describe the high-resolution radio observations and subse- 
quent morphologically classification of the sources in the sample, 
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along with the subsequent space density measurements of the de- 
tected FRIs. 



ACKNOWLEDGEMENTS 

EER acknowledges a research studentship from the UK Particle 
Physics and Astronomy Research Council. PNB would like to 
thank the Royal Society for generous financial support through 
its University Research Fellowship scheme. The United Kingdom 
Infra-red Telescope is operated by the Joint Astronomy Centre on 
behalf of the UK Particle Physics and Astronomy Research Coun- 
cil. We acknowledge the UKIRT Service Programme for some of 
the infra-red imaging. The Isaac Newton Telescope is operated by 
the Isaac Newton Group, and the Telescopio Nazionale Galileo is 
operated by the Centra Galileo Galilei of the Consorzio Nazionale 
per l'Astronomia e l'Astrofisica, both in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astrofisica de Ca- 
narias. The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities Inc. 



REFERENCES 

Barbara G. & Poggianti B. M, 1997, A&A, 324, 490 
Baum S.A., Zirbel L., O'Dea C. P., 1995, ApJ, 451, 88 
Bershady M.A., Hereld M., Kron R.G., Koo D.C., Munn J.A., Ma- 
jewski S.R., 1994, AJ, 108, 870 

Best P. N, van Dokkum P. G, Franx M., Rottgering H. J. A., 2002, 
MNRAS, 330, 17 

Best P. N., Kauffmann G, Heckman T. M., Ivezic Z., 2005, MN- 
RAS, 362, 9 

Best P. N, Kaiser C. R., Heckman T. M., Kauffmann G, 2006, 
|astro-ph/0602171| 

Blanton M. R. et al., 2001, AJ, 121, 2358 

Bower R. G, Benson A. J., Malbon R., Helly J. C, Frenk C. S., 

Baugh C. M„ Cole S., Lacey C. G, 2006, MNRAS, 370, 645 

Bruzual G, Chariot S., 2003, MNRAS, 344, 1000 

Casali M. & Hawarden T, 1992, UKIRT Newsletter, 4, 33 

Clewley L. & Jarvis M., 2004, MNRAS, 352, 909 

Condon J. J. & Yin Q. F, 1990, ApJ, 357, 97 

Dunlop J.S. & Peacock J.A., 1990, MNRAS, 247, 19 

Eales S., Rawlings S., Law-Green D., Cotter G, Lacy M., 1997, 

MNRAS, 291,593 

Fanaroff B.L and Riley J.M., 1974, MNRAS, 167, 31 
Gawronski M. P., Marecki A., Kunert-Bajraszewska M., Kus A. 
J., 2006, A&A, 447, 63 

Gopal-Krishna & Wiita P.J., 2000, A&A, 363, 507 
Jackson C.A. & Wall J.V., 1999, MNRAS, 304, 160 
Jorgensen I., 1994, PASP, 106, 967 

Kauffmann G, Heckman T. M., Tremonti C, Brinchmann J., 

Chariot S., White S. D. M., Ridgway S. E., Brinkmann J., 

Fukugita M., Hall P. B., Ivezic Z., Richards G. T, Schneider, D. 

P. 2003, MNRAS, 346, 1055 

Koo, D. C. & Kron, R. G, 1982, A&A, 105, 107 

Kron, R. G. 1980, ApJS, 43, 305 

Landolt A. U, 1992, AJ, 104, 340 

Ledlow M.J. & Owen F.N., 1996, AJ, 112, 9 

Leggett, S. 2005, http://www.jach.hawaii.edu/UKIRT/astronomy/ 

calib/phot_cal/cam_zp.html 



McCarthy P. J., van Breugel W., Spinard H, Djorgovski S., 1987, 
ApJ, 321, L29 

Monet, D. G. et al. 2003, AJ, 125, 984 
Oort M. J. A. & van Langevelde H. J„ 1987, A&AS, 71, 25 
Oort M. J. A., Katgert P., Steeman F. W. M., Windhorst R. A., 
1987, A&A, 179, 41 

Oort M. J. A. & Windhorst R. A., 1985, A&A, 145, 405 

Owen F.N., Laing R.A., 1989, MNRAS, 238, 357 

Perley R.A., 1999, in Synthesis Imaging in Radio Astronomy II, 

ASP Conf Ser 1 80, 383 

Taylor M., 2000, http://www.star.bris.ac.uk/~mbt/wfc/calibrate/ 

Sandage A., 1972, ApJ, 173, 485 

Smith J.A. et al., 2002, AJ, 123, 2121 

Snellen I. A. G. & Best P. N, 2001, MNRAS., 328, 897 

Snellen I. A. G, Bremer M. N, Schilizzi R. T, Miley G. K, van 

Ojik R., 1996, MNRAS, 279, 1294 

Stoughton C. et al., 2002, AJ, 123, 485 

Wall J.V., 1980, Phil. Trans. R. Soc, A296, 367 

Waddington I., Windhorst R.A., Dunlop J.S., Koo D.C., Peacock 

J.A., 2000, MNRAS, 317, 801 

Waddington I., Dunlop J.S., Peacock J.A., Windhorst R.A., 2001, 
MNRAS, 328, 882 

Willott C. J., Rawlings, S., Blundell K. M„ Lacy M., 1999, MN- 
RAS, 309, 1017 

Willott C.J., Rawlings S., Blundell KM., Lacy M., Eales S.A., 
2001, MNRAS., 322, 536 

Willott C. J., Rawlings, S., Jarvis M. J., Blundell K. M„ 2003, 
MNRAS, 339, 173 

Windhorst R.A., van Heerde G.M., Katgert P., 1984, A&AS, 58, 
1 

Windhorst R. A., Mathis D., Neuschaefer L., 1990, in R. G. Kron 
(ed.), ASP Conf. Ser. 10: Evolution of the Universe of Galaxies, 
pp389-403 

Windhorst R. A., Fomalont E. B., Partridge R. B., Lowenthal J. 

D„ 1993, ApJ, 405, 498 

York D.G et al., 2000, AJ, 120, 1579 

Zirbel E.L. & Baum S.A., 1995, ApJ, 448, 521 

This paper has been typeset from a TgXy KTgX file prepared by the 
author. 



APPENDIX A: IMAGING AND SPECTROSCOPY NOTES 
ON INDIVIDUAL SOURCES 

Al Imaging source notes 

53w054a and 53w054b - Windhorst et al. (1984) classified these 
two sources as radio galaxy lobes but, in agreement with Wadding- 
ton et al. (2000), optical identifications were found for both sources 
indicating they are two separate radio galaxies. However, the 
Waddington et al. identification of 53w054b seems incorrect - a 
faint source is detected in the A"-band which is more closely asso- 
ciated with the radio position. 

53w087 and 55wl32 - Neither of these sources were detected in 
the optical or infra-red imaging and, as a result, are quoted as be- 
ing at z >4 in Table [14] However, inspection of their radio maps 
shows that both appear to be extended, diffuse-like sources, which 
suggests that they are more likely to be dusty star-bursting galax- 
ies, or post-mergers, at z = 1-4. Their estimated redshift limits 
may, therefore, be unreliable. 
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66w042 - Whilst the centre of this source does not appear to align 
with the indicated optical galaxy, the identification is valid as the 
faint radio core is offset to the west. 

55wl33, and 55wl43a/b - The r-band magnitudes measured for 
these 3 sources are all around the la level. They therefore should 
be treated as unreliable. 

55wl37 and 60w032 - The slightly different telescope pointing 
used in the 2004 INT observations, compared to that used in 2003, 
meant that these sources were not present on the 2004 images. 
55wl59b - The magnitude measured for this source should be 
treated as less reliable due to the presence of two other objects in 
close proximity to it. 

A2 Spectroscopy source notes 

53w070 - The single broad line is identified as Mgll, giving a red- 
shift of 1.32. 

53w086b - The 4000A break can be seen in this spectrum but no 
lines are convincingly detected. The redshift of 0.73 is consistent 
with the r-band magnitude estimate of 0.81. 
53w089 - Waddington et al. (2000) detect [Oil] and [OIII] for this 
source, giving a redshift of 0.635 but nothing is detected in these 
(shallower) observations. 

66w031 - The identification of single strong line as [Oil] at a red- 
shift of 0.81 is consistent with the redshift estimated from the r— 
band magnitude (0.97). 

66w036 - No line data are quoted for the G-band absorption line 
for this source as the negative extent of the flux density suggests 
that it is contaminated by noise. 

55wl24 - The single broad line is identified as Mgll, giving a red- 
shift of 1.34. 

55wl28 - The weakly detected single line is identified as [Oil], 
giving a redshift of 1.19, but this is uncertain. 
55wl31 - The single line is identified as [Oil] giving a redshift of 
1.12, which is consistent with the redshift of 1.15 estimated using 
the r-z relation. 

55wl40b - The redshift of 1.68 is from a single, broad line, identi- 
fied as Mgll. 

60w039 - The two observations of this source on masks L3 and L4 
have been combined together with the scomb task in iraf. 

APPENDIX B: THE RADIO/OPTICAL IMAGES 
APPENDIX C: THE SPECTRA 

Here Figure [Cl] shows the spectra resulting from the DOLORES 
MOS observations. 



Figure Bl. See attached jpg files. The identifications and VLA A-array 
contour maps for the Hercules field; from left 1 .4 GHz radio contour map, 
r with overlay if present, i and K with overlay if there is no r detection. 
Where necessary, host galaxy positions are marked with crosshairs and 
some images have been gaussian smoothed for clarity. Radio contours start 
at 24/xJy/beam and are separated by factors of \/2. The primary beam cor- 
rection has not been applied to the radio maps so that uniform images can 
be presented. 

Figure B2. See attached .jpg files. Similarly, the identifications for the Lynx 
field. From left, 1 .4 GHz radio contour map, r with overlay if present, i and 
K with overlay if there is no r detection. Radio contours start at 24/iJy/beam 
and are separated by factors of \/2. The primary beam correction has not 
been applied to the radio maps so that uniform images can be presented. 

Figure B3. See attached .jpg file. The un-primary beam corrected radio 
contour maps for the sources in both fields, in the complete sample, with no 
optical or infra-red identification. Radio contours stall at 24/xJy/beam and 
are separated by factors of 

Figure B4. See attached jpg file. The radio images and infra-red identi- 
fications, if present, for the sources not included in the complete sample. 
Radio contours start at 24/iJy/beam and are separated by factors of V2. 
The primary beam correction has not been applied to the radio maps so that 
uniform images can be presented. 
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Figure CI. The spectra resulting from the MOS observations. Residual sky features are marked with a and dashed lines indicate 
the position of the 4()0()Abreak. 2D spectra are also shown in cases where the line detection is weak. 
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Figure CI - continued 
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Figure CI - continued 
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Figure CI - continued 
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Figure CI - continued 
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Figure CI - continued 
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